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Different possible domain structures of single crystals of nickel and iron have 
been considered. In each case the crystal has the shape of a rectangular rod with the long 
axis in the [100], [110] or [111] direction. The influence of an external magnetic field 
directed parallel to the long rod axis has been examined. The structures found consist of 
main domains and surface domains closing the magnetic flux. Domain structures of 
nickel crystals are found to be analogous to those of iron. 


1. Introduction 


The aim of this paper is to investigate the domain structure of some nickel and 
iron single crystals. The mode of reasoning is essentially the same as in the known 
papers of Landau and Lifshitz (1935), Néel (1944) and Lifshitz (1944), who investigated 
the domain structure of cobalt and iron crystals. For single crystals of nickel, no 
calculations of this type have been published. We shall'also briefly discuss the in- 
fluence exerted on the domain structure by a.demagnetizing field near the ends of 


a crystal rod. : 


2. Method of Calculation 


We assume the crystal to have the shape of a restangular rod: of infinite length 
and with its axis in [100], [110] or [111] direction. The surfaces of the crystal are 
cut symmetrically with respect to the crystal cubic axes. One pair of crystal planes is 
assumed large compared with the other: crystal width L2» crystal height h. We then 
have six possible cases for nickel and six for iron crystals, some of them being equi- 
valent. ! d : 
The probable domain structure can be guessed on the basis of symmetry canisi- 
derations. This procedure is simplified when a uniform external magnetic field is 
assumed to act, since it favours some of the easy directions in the crystal. We then 
write down the different contributions to the free energy of the crystal and minimize 
them with respect to some parameters which determine the size of the structure. 


(169) 
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3. Contributions to the Domain Energy 


We shall use the following energy density expressions: 
1. Magnetocrystalline energy: 
f, = Ky(aiad + ajaz + a3c1) (3.1) 


where K, is the anisotropy constant, and œ; are the direction cosines of the magneti- 
zation vector in a domain. 


2. Magneto-elastic energy: 
Sme = — 368 Ayoo (af €i + 063829 + 063835) — 3C3 Ain * 2 (06955615 + 6591653 + 6505831) + 
+ $ Calen + E22 + E38)? + Co(et + 622 + 633) + 2Cs(ei2 + €23 +31) (8-2) 
where A499. 4,1 are magnetostricion coefficients, &,, are components of the strain 
tensor and the constants C,, Ca, C4 are given in terms of the elastic coefficients c, 
ciz» Can» by: 
Ci = Cy, Ca = $ (eu — Cy), Cy = Cu. 


3. Bloch wall energy: 
9, 
Ja = eves © TU do (3.3) 
9, 


where y, is some constant depending on the exchange energy for parallel spins, the 
polar angles © and ® determine the direction of the magnetization vector I in the 
wall and f, denotes the excess magnetocrystalline and magneto-elastic energies in 
the wal' divided by the modulus of the anisotropy constant K, (Lilley 1950). 


4. Self energy of a permanent magnet in its own field: 


> -= 


Um PEU. il: H, (3.4) 
where H, is the magnetic field in the domain, and J the magnetization vector. 


5. Energy of interaction of a permanent magnet with a uniform external magne- 
tic field: 
T DIT (3.5) 


with H as external field. 


4. Domain Structures of Nickel Crystals 


A. Crystal rod long axis [100], main plane (011). An external magnetic field H 
in [100] direction, privileges the [111] and [111] easy directions. Therefore the magne- 
tization vectors in the main domains will lie in these directions when the field is small, 
and will move towards the field direction with rising field intensity. The expected 
domain structure with the closure domains removing the magnetic poles is shown 
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in fig. (la, b). Structure (a) has been suggested by Bates and Wilson (1953), structure (b) 
differs by another shape of the Q domains. Both structures resemble that of Néel 
(1944) for an iron erystal, and the following analysis is analogous to Néel's. 

The double domain width d which determines the size of the structures is obtained 
by minimizing the sum of the excess magnetocrystaline, magneto-elastic, and external 
magnetic field energies in the closure domains and energy of the wall between main 


Fig. la. Domain structure (a) of a nickel crystal when the long axis is parallel to the [100] direction. 
The magnetization vector in the Q domains is parallel to the [110] direction. 


Fig. 1b. Domain structure (b) of a nickel crystal as in the preceding case. The magnetization vector 


in the Q domains is parallel to the pij or nl direction. 


domains. The energy of triangular areas of magnetic poles on the (011) surfaces and 
of the walls of the closure domains is small compared with the other terms and may 


be neglected. 
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Domain Structure (a) 
Using (3.1) and (3.5) we shall calculate the expressions for excess magnetocrystal- 
line and external magnetic field energies in the Q and P closure domains (fig. 2). 


We have 


AfL = 1 K, (3 cos?9 — 2 cos?O) (4.1) 
Afg = — K, (3 cos*@ — cos*@) (4.2) 
Af? = — 4 Ky (1 +2 cost — 3 cost6) et 
Afi = — K, (1 — 3 cos? O) (cos © — cos? O) (4.4) 


The external field H has been eliminated with help of the equation 
= z (cos 9 — 3 cos? O) 


obtained from the condition of minimum magnetocrystalline and external magnetic 
field energy of the crystal (Stewart 1955). 

In evaluating the magneto-elastic energy term we shall follow Kittel (1956) and 
give only its upper limit. The closure domains form triangular prisms. Since magneto- 


(1-x) a 


Fig. 2. Shape of closure domains of structure (a). The magnetization vector in the Q domains is perpendi- 
cular to the plane of the paper. 


striction causes shortening of the length of every domain, domains magnetized in 
different directions will not fit together smoothly except by expenditure of elastic 
energy in forcing the material to fit together. Thus, the domains of closure may be 
regarded as being elongated to join on to the main domains. The strain imposed by 
the main domain structure is given by the equation 


3 l 
7-5 rm (a Pt + a3 Bà + as B3 — = + 3Ayı (a 25 By b2 + 
+ az Gs Bz Bs + as ay Bs B1) 
where a; are the direction cosines of the magnetization vector and f; the direction 
cosines of Öl. 
For Q domains the strain in the [011] direction is equal to: 


él 
DA MARO 
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and the significant components of the strain tensor are: 


ôl 


£99 = Ey = 7 


the remaining ones vanish. Thus for the excess elastic energy ofa Q domain (the energy 
of basic structure being zero in this approximation) we have 


t 
Af? = Te (en + c12) (A100 + 3443)? (4.5) 
For P domains the strain in [100] direction equals to 


ôl 


TPP ^10 = Eu 

the remaining components of ej, vanish. The excess elastic energy term is now: 
mou 2 

Afa = 3 ĉn A100 (4.6) 


The shapes of Q and P domains are determined by the condition that the normal 
component I, of magnetization vector be constant across the wall (fig. 2). The corres- 
ponding volumes per unit of crystal height are 


vg = 2 d? tan ©, üp = ! a — x)? d? cot "d (4.7) 


RI 
eI 


We shall now write down the energy equation. There are 2/d closure domains 
of every type (per unit length) on both sides of the crystal and 2/d walls with energy fọ 
per unit area, on unit length of the crystal. The sum of the excess energies and the 


wall energy is given by: 


Fa E [33 d? Wo + (1 — 3)? d? Wp + [fal (4.8) 
where | 
2008 
Wo = s tan (AR + Af + Af, = aq ets (ASE Aff Af) 


and L is the crystal width. l 
Minimizing F with respect to x and d we obtain for the double domain width d 


] / fa L (Wo + Fp) 
a= 7v: (4.9) 


which has the form of Néel's formula for the analogous structure of an iron 


the equation: 


crystal. 
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Domain Structure (b). 
This differs only by Q domains. The corresponding excess energies now have 


the following form: 


1 1 
Aft =7 [s — 2 cos? O 306] (4.10) 
1 
Af = K, cos O bs + cos o) (1 — 3 cos? ©) (4.11) 
1 
Af = Zh (ei. + 63) Ain (4.12) 


As regards the last expression the Q domains are assumed to be elongated in the 
[011] direction with strain $4, and strain tensor components Es = £33 = $ Agi, 


xd (1-x) d 


Fig. 3. Shape of closure domain of structure (b). The angle between the magnetization vector in the Q 
domains and the plane of the paper is « = arc cos 1/Y3. 


the remaining ones being equal to zero. For the domain volume per unit heigth we 
obtain the expression (fig. 3) 


1 sin © 1 
UO — ad rp oer E 
4 cosO +1jy3 4 


As in preceding case double domain width d is given by the formula: 


1 / fa L (Wo + Wy) 
d= P = P (4.14) 


with sin O/(cos © + 3-4) instead of tan O in IV g: Numerical results for both struc- 
tures (a) and (b) are collected in tab. 1, and plotted in fig.4. Experimental points are 
taken from the paper of Bates and Wilson (1953). The formula for wall pigra Fy is 
calculated at the end of the paper. 

Structure (b) is associated with a lower energy and one may m find this 
structure rather than structure (a). 

The quantitative disagreement with experiment is due mainly to the neglect of 
the demagnetizing field near the ends of the crystal, as has been shown for an iron 
crystal of analogous structure in the author's previous paper (1958). For the main 
plane (011) the same structure results. 


x? d? tan 9 (4.13) 
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Tab. 1 
1 z 
cos 0 = I/I, | H (Oe) V6 (erg. cm~?) W d(erg.cm-) Wp(erg. cm”? |f, (erg.cm"?)| d, (cm) d; (cm) 
1/3 0 2.56 x 103 0.07 x 103 9.65 x 10? 0.17 921025 DEDI S 
0.64 16 3.50 1.45 9.03 0.14. 7.5 10.6 l 
0.67 24 4.00 21 9.94 0.12 6.6 8.4 
0.70 34 4.55 2.86 8.56 0.11 6.1 7.2 
0.76 58 5.62 = AR 7.72 0.09 5:3 5 
0.81 82 6.40 5.47 6.86 0.06 4.3 4.4 
0.90 133 6.19 6.10 5.33 0.02 2.6 Dell 
0.94 159 6.14 6.36 5.55 0.01 1.9 1.8 


Tab. 1. Values of d for structures (a) and (b) calculated by means of formula (4.9) and (4.14), with values. 
of the constants as given in last sections of the paper. 


Crystal rod long axis [110], main plane (110). 
In an external magnetic field parallel to crystal rod axis the domain structure is 
expected to be as shown in fig. 5. In a weak field the main domains are magnetized 


in the [111] and [111] directions, and with increasing field strength the magnetization 
vectors rotate towards the field direction. The magnetic flux is completely closed by 
triangular closure domains when the angle v satisfies the condition tan y — cot (0/2) 


where O is the angle between the magnetization vector / in the main dómains and 
field direction. 


20 
16 


72 


120 160 


40 80 
H oersteds | 
i in wi ith i i ] magnetic field H. The curves cor- 
ig. 4. Ch f double domain width d with increasing external } r 
> SPEER with L = 0:8 cm as in the paper of Bates and Wilson (1953), from which the experi- 
1 mental points have been taken. 
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Fig. 5. Domain structure of a nickel crystal with the long axis parallel to the [110] direction. 


The minimum condition is formed by means of excess magnetocrystalline, elastic, 
and external field energy in closure domains and wall energy between main domains. 
These contributions have the following form: 


Aff Ky : — cos? O + 7 cos? o) (4.15) 

Afu = K, cos Ø (cos 9 — 1) (2 — 3 cos? O) (4.16) 
5 1 

Afa = 39 (ci + C12) (99 + 343)? (4.17) 


So = 2y, sin O | sin O i — 2 sin? ej -+ 


+ (V3/2 — /2/3 - sin? ©) arc sin DONEC edt (4.18) 


1 2 
sn?O 3 


The last expression is calculated at the end of the paper. The volume of a closure 
domain per unit crystal height is 4 d? cot (0/2). 
For double domain width d we obtain 


d = Vf, LIW (4.19) 


with W = 2 cot (0/2) (Af, + Afy + A If). In tab. (2) numerical results are collected. 
The corresponding is plotted in fig. 6. 

For the main plane (001) the same structure results. 

For short crystals the influence of the demagnetizing field should result in in- 
creasing the double domain width d as in the preceding case. 
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Tab. 2 


W (erg. cm) 


cos © = IJI, N H (Oe) Ju (erg. em?) 


¥2/3 5.65 x 103 : 
0.86 5.14 5.5 
0.90 4.93 4.8 
0.92 4.93 4.4 
0.94 5.17 3.8 
0.95 5.40 36 


Tab. 2. Values of d calculated by means of formula (4.19) with values of the constants as in the 
preceding case. 


Crystal long axis [111], main plane (112). 

The crystal favours the easy directions lying in the (110) plane so that the expected 
domain structure should be as in fig. 7. A weak magnetic field in the [111] direction 
should magnetize the crystal to saturation. The magnetic flux is completly closed with 
a + B = 2/2 and sin a = 1/Y3. The domain width of the equilibrium configuration 
is determined by the elastic energy of the closure domains and the wall energy of the 
main domains. The volume of closure domains per unit area is given by dy2/6. To 
find the elastic energy term we assume an elongation of closure domains in the [111] 
direction forced by the main domains. The resulting elastic energy density is equals to 


D s 
Af. "pe 3 tu Ann (4.20) 


60T a x 105 (em) 


5.0 


4.0 


3.0 60 


20 40 
H(0e) ; ps. > 
Fig. 6. Change of double domain width d with increasing external magnetic field H, according to tab. 2. 
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The wall energy per unit area of the (112) surface amounts to 


oe : ar Im: (4.21) 


where h is crystal height. For domain width d we find 


= 5 Vifalfa (4.22) 


With h = lcm, f, = 0.32 erg cm? taken from the paper of Lilley (1950), we 
obtain d = 2:3 X l0 4em. 


Fig. 7. Domain structure of a nickel crystal with the long axis parallel to the [111] direction. 


5. Domain Structures of Iron Crystals 


For the crystal rod long axis in the [100] direction, the domain structure found by 
Lifshitz (1944) exists. The two planes (010) and (001) are equivalent. 

When the crystal long axis points in the [110] direction and the main plane is the 
(001) one, we have Néel’s domain structure and for the (011) main plane, the domain 
structure discussed by Bates and Hart (1956). 

We will now find the domain structures of crystals with the long axis in the [111] 
direction. 

Long axis [111], main plane (110). 

In a very weak magnetic field in the [111] direction, the magnetization vectors 
of the main domains are distributed among the [100], [010], and [001] easy directions, 
and they approach the field direction with growing field intensity. This is in accor- 
dance with Parker’s analysis of iron single crystals by an electrical method (Parker 
1956). The magnetic flux on the main crystal surface is completely closed, the angle v 


between the inclined part of the wall and the (112) surface being determined by the 
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condition tan y = tan V3/y2 fig. 8. On the (112) crystal surface some domains of 


closure should appear to diminish the magnetic energy, but their structure is difficult 
to predict. 


The size ofthe structure determined by the width of two domains d, +d, — d 


is obtained by minimizing the magnetic energy of the (112) surfaces and the wall 
energy. The corresponding energies per unit area have the form 


9 JE d sin? @ 
ee ear XE — sin? e e d (5.1) 


(5.2) 


Fig. 8. Domain structure of a iron crystal with the long axis parallel to the [111] direction and. the main 


plane (110). The value of d, is always 2d, which is not represented by the model. 


Tab. 3 


cos 0 = IJI, 


3.23 x 10° 1.52 x 10-3 


0.67 2.66 1.53 
0.76 2.04 1.53 
0.81 1.66 151 
0.90 0.92 1.44 
0.94 0.56 1.34 


Tab. 3. Values of d calculated by means of formula (5.3) with values of the constants as given in last 
sections of the paper. 


These formulae are calculated in the last part of ihe paper. 
The minimum condition p 


; = (Lf, | Mf E (5.3) 
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where L is the crystal width. Numerical results are tabulated in tab. (3) and plotted 
in fig. 9. 

The experimental structure should be characterized by a larger value of d because 
of the reduction of M by means of closure domains. 

For short crystals the influence of the demagnetizing field will be important. 
This is connected chiefly with the unknown closure domains and therefore cannot be 


predicted. 
Crystal long axis [111], main plane (112). 
d x 105 cm 
15 
1.2 
19 
0 50 100 150 200 250 
H (Oersteds) 


Fig. 9. Change of double domain width d with increasing external magnetic field H, according to tab. 3. 


The main domain structure is the same as in the preceding case; but now the 


(112) surface is more important and the magnetic flux is closed by means of triangular 
closure domains magnetized in the [111] direction, fig. 10. It is not possible to remove 


simultaneously the magnetic poles from the (110) surfaces. The closure domain struc- 


Fig. 10. Domain structure of an iron crystal with the long axis parallel to the [111] direction and the main 
plane (112). 


ture which exists actually will lower the magnetic energy of the configuration. The 
calculated energy of alternating strips of magnetic poles presents only an upper limit, 
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The width d of two domains is determined by the minimum condition for the 
following energy terms: 


The excess magnetocrystalline energy in closure domains 


1 
di eK, E — sin? (O, — O) + Ž sin’ (9, — e) (5.4) 
- The magnetic energy of the strips of poles (calculated at the end of the paper): 


41? d l| 2 1 
E = E (9, — Oe Và sin (O, — e| . Ds E err )- Md (5.5) 
odd integers 


_ The excess external field energy in closure domains: 


l 
| pa a ur |n (20, — 20) (a - È sin? (©, - eJ] (5.6) 


_ The Bloch wall energy of walls between main domains: 


"de = zig iE n 20)? (5.7) 
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For the width d of two domains we obtain: 
f= | 2h jé (5.8) 
MIL 3- S(Af. + Afi] 
where L and h are the width and height of the crystal. 
Numerical results are tabulated in tab. (4) and plotted in fig. 12. 


Tab. 4 
cos O = IJI, |H(Oe)| M (erg. cm?) iy Afp (erg. cm 3)| f, (erg. em?) d (cm) 
1/y3 0 1.68 x 105 0.93 x 105 0.73 91X I0 
0.67 55 6.35 0.61 0.62 9.2 
0.76 118 4.86 0.32 0.48 9.4 
0.81 156 3.98 0.20 0.38 9.3 
0.90 224 2.20 0.04 0.19 9.1 
0.94 250 1.34 0.01 0.10 8.5 


Tab. 4. Values of d calculatedby means of formula (5.8) with values fo the constants as in the preceding case. 


d x 10 Som 


1.0 


0.8 


0.6 
0 50 700 750 200 250 
H (Oersteds) 


Fig. 12.-Change of double domain width d with increasing external magnetic field H, according to tab. 4. 


The experimental values of d should be larger since we have overestimated the 
magnetic energy term in the denominator and because of the demagnetizing field 
influence. 


6. Calculation of Bloch Wall and Magnetic Energies 


Bloch wall energies for nickel crystals. 
We will use the expression (3.3) for the wall energy per unit area 


9, 
fa = 2ygsin 8 | fido | (6.1) 
9, 


For the [100], (110) crystal, the wall lies in the (100) plane. The magnetization 
vector ve in the wall, turns from direction the (©, ®,) to the (©, 9) (fig. 13). The 
angle @ is constant across the wall, is equal to ©, = arc cos 1/y3 in zero external 
magnetic field, and diminishes with growing field strength, while ® varies from 0 
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to x. Expressing the direction cosines of L by 


© and © we obtain for the excess magne- 
tocrystalline energy the expression 


wa ain al i 
f: u — fim = Af, = 44 sin* @ sin? 26 (6.2) 


In evaluating the magneto-elastic energy term we shall follow Lifshitz (1944). 
We need the values of the strain tensor £j. In the domains, for x = + oo (fig. 13) 


the magnetization is uniform and magnetostriction is also uniform; therefore all the 
components of the stress tensor c; are equal to zero. Using the relation 
= Of me 

Or 
we obtain the components of the strain tansor £j, in the domains. The direction cosines 


= 
of I, are 


Oik 


(6.3) 


l 1 
æy = cos O, a= t ygn O, Ge EE E ee 


and the components of ¢, in the domains are 
given by: 


o 


[3 


C 
B Aes (cos O ar) 


3 1 C 
== Seren — sin? ET ee a 
£j = Es = 5 A00 E sin? @ 30 42 a 
3 PN 
fgg = — q ^u sin? 9, 


3y2 A 
£g = — Eg = T 


Since the distribution of magnetization in 
the wall depends only on x coordinate, the strain 
tensor £j, should also depend only on x. It may fig. 13, The angles O and Ø in the (100) 
then be deduced from equations defining the Bloch wall of a nickel crystal. 
components of strain tensor £j, that £55, £53; £33 ‘ 
must be independent of x and equal to those in the domains. Further, the 
equilibrium condition of an elastic body is given by 

905 
9x, 


== 0 


Since the components of the stress tensor o; depend also only on x, this yields the 


relations: 
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But oj are 0 for x = + © and therefore also in the wall 
Opx = 0,, = Oy = O 
Condition (6.3) now determines the remaining components of the strain tensor: 


C 


3 3 . o 
en = y A100 (cos 0 — sc gn] | =T Ai sin 20 - cos (45° — ®) 


E13 = TW" T sin 2 0 * COS (45° + 9) 


where the relations: 


a, = cos O, a,=sin@-cos(45°—), a, — — sin O - cos (45° + o) 


z 


for the direction cosines of J, in the wall have been used. We may now write the 
expression for the excess magneto-elastic energy: 


nir D f fomen p A 4244 sin‘ © - sin? Ø (6.4) 
The term f, = [Kil (Af, + 4/,,,) is equal to: 
fa = sint © sin? Ø (cos o — zit | A) - (6.5) 


After inserting (6.5) in (6.1) we obtain for the wall energy: 


2 
(1 + a)? —1 
with a = 9C4A2,/2 |K,]. We take for yọ the value given by Lilley (1950), yo = 
= O.l4erg- cmr?, the other constants are: C3 = cy = 1.185 x 10!2 erg cm? 


fy = 2 Yq sin? O |a +a -- aln 


(Bozorth et al. 1949), K, = — 0.5 x 101? erg cm”? (Stewart 1955), 444 = — 19 x 107$ 
(Corner et al. 1955), Thus finally 
So = 0.3F sin? O (6.6) 


which has been used in tab. (1). 

For the [110], (001) crystal the wall is in the (110) plane, fig. 14. The magnetiza- 
tion vector 1, turns from the (O, O) to the (©, x) direction. In a zero external magnetic 
field © = ©, = arc cos /2/3. The excess magnetocrystalline energy is in this case 


equal to: 
foal domain _ Af, = — K, sin? © sin? Ø atapse - 3 sin) 
(6.7) 


In order to calculate the magneto-elastic energy term we pass to the x’y'z’ system of 
axes (fig. 14). Now the wall normal is parallel to x’ and we may repeat the foregoing 
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considerations. In the x'y'z' coordinate s 


ystem the magneto-elastic energy expres- 
sion has the following form: ae 


M l 1 
Sime = — 3 A, 99 EE (€11 — 2&2 + E22) + p &3 (E11 + 26g + E22) + x? a E 
SGAM. 2 : M : 

3 ^11 3 Oy & (E11 — E22) + y2 5 Xz (£i; + E93) + y2 Hs 04 (E13 — en) "E 


1 1 1 
HE p C, (En + €g9 + Eza)? + Cy E (£u + 63)? + 281, ai a + 


il 
+ 26s E (£i — £29)" + e2, + a |. 


The a, and £ refer to x'y'z' axes. 


Fig 14. The angles © and @ in the (110) Fig. 15. The angles O and Ó in the (111) 
Bloch wall of a nickel crystal. Bloch wall of iron crystal. 


After having found the expressions for ej, in the domains and in the wall we obtain 
for the excess magneto-elastic energy the formula: 


DNI as, 
Cy + Cg + 63 
+ egs [ab + 3 sin?@ (e, + ca + ¢3) C2 A19] + 


p 
as E sim? 20 - (c, + cg + c) (ca Afi — % A) ee DEM UCR sche TS 


A sin? Ø + B sint Ø 6.8) 
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where £54, & 3 are components of strain tensor in the domains and 
2 Ze 
red oy cos? © (C; Ayoo + Ca An)» pe 3 sin? © (C, Ayoo — Cs À111). 


After inserting the numerical values for the constants we obtain for 4 and B: 
A = sin?0 - (2.5 cos?@ —24.0) x 10?, B= —sin* 0- 2.85 x 10? 


Comparison with coefficients of sin? and sint Ø in (6.6) shows that A and B 
may be neglected in calculating the wall energy. The values of the constants used are:: 
Aioo — —52 x 10-9 (Corner et al. 1955), C= $(e—05) with cy = 2.50x 
x 1012 ergcm”?, C = 1.60 x 10! erg cm è, C, = cy = 1.185 x 1012 erg cm? 
(Bozorth et al. 1949). 

With help of (6.1) and (6.6) we find the formula for wall energy: 


To = 2y, sin O | sin © (3 — 2 ie]? + 


1 
1 
sin? © 
which has been used in calculating the values of tab. (2). 
Bloch wall energy for iron crystals. 
For both of the crystals considered the wall lies in the (111) plane and the magneti- 
zation vector i turns from the (©, O) to (©, 27/3) direction (fig. 15). In zero external 


magnetic field © = Oy = arc cos 1/|/3. In the case of iron the magneto-elastic energy 
may be neglected. The excess magnetocrystalline energy Af, is given by 


+ (73/2 — y2]3 ` sin? ©) arc sin (6.9) 


Bu 
3 


Af, =K, e sin? O cos O (1 — cos 3 Ø) (6.10) 


By means of (6.10) we obtain for Bloch wall energy the expression 
fa= $ yosin? O e Ee eJ (6.11) 


which has been used in tab. (3.4). For yọ the value 0.62 erg cm-? has been taken 
(Lilley 1950), and K, = 2.8 x 105 erg. cm”? (Williams et al. 1949). 


Magnetic energies 


We shall consider the magnetic field energy of coplanar strips. of poles following 
Kittel (1949). Let the strips lie in the xy plane, the y axis being parallel to the axis 
of the strips. In order to find the magnetic field H of the two-dimensional problem 
we shall expand the surface density ø of magnetic poles in a Fourier series Yo, and 
then write down the solution for potential ® in the form © (x, z) = 2X4A,o, exp 
(—nz|z||l) with A, determined by the boundary condition. 
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For the [111], (110) crystal, the surface density ø is expressed by a symmetric 
function (fig. 16) with the period d = 2 (dj + dj). When expanded into a Fourier 


series it is equal to 


2 STET eov 
g — 3 (hd, — I, dj) 4 er aan (2m $) cos (2an &) 


d 


Fig. 16. Surface density of the magnetic poles on the (112) plane represented by a periodic function. 


and the potential ® has the form 
EEUU ex 1 d, x Iz] 
P= A PISO 2 74 A, sin | 2zn g] 2zn j] xP \ 2zn F (6.12) 


The boundary condition [9D/dz]_» —[9G/9z] 45 = 4mo yields: 
v res s I d, = 1, d,. 
n 
The last condition determines the relation d,/d, for the domain widths. With the help 


of (6.12) the magnetic field energy per unit area —1/2dl- i I- Hdv (l is the length 
of a strip) is given by 


rE d,\ E R E E 1. N „. 
cope EE n 2an 7 er Gia, pase 3] 619 
1 1 j 


mag — e" Agr? 


after inserting J, = I, sin O, I, = $/,sin@, from which d, = 2d, follows. The 
approximate value of the sum in (6.13) is 1.14 which is smaller than the limit value 
with error < 0,005. Thus finally x 


Snag = 0485 x 10° d sin? O 


In the case of [111], (112) crystal, the surface density o is given by a symmetric 
function of period 2d = 2 (2d, + d,) (fig. 17). For potential ® we obtain the expres- 


sion 


Ø = Al = A, sin (am) cos Ea cos (>=) exp (—nz |z|jd) (6.14) 
1 , 
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with A, = 2d/n obtained from the boundary condition, The magnetic energy per 
unit area is equal to 


(6.15) 


8I? d loe? es 


— COS 
mag = - n3 ad 


odd integers 


After inserting I = I2 > [cos (0, — ©) — 1//2- sin (Oy — ©)] and d,/2d = 
= tan y,/2 (tan y, + tan y,) = 1/6 (fig. 11), we obtain 


f(x) 


4 
JEE = 


Id * 1 
[cos (0, — 9) — 1/2 - sin (9, — 6)]? - 5 z cos? (=) 
odd integers 
with ©, = arc cos 1//3. To obtain the approximate value of the sum we take the 
first six terms with the error < 0.005. 


Thus finally: 


2 
Scag = 0.168 E (0, — 9) — a sin (0, — e| d x 10° 


yà 


which was used in calculating tab. 4. 


Conclusions 


The domain structures of nickel crystals are found to be analogous to those of 
iron crystals. The domain structure of the [111] nickel crystal corresponds to that 
of the [100] iron crystal, the domain structure of the [100] nickel crystal to that of 
the [110] iron crystal and domain structure of the [110] nickel crystal to the [111] iron 
crystal. The external magnetic field influences the domain structures in similar way 
causing a decrease of domain width d, except for the [111] iron structures where the 
value of d at first increases with increasing field strength. 


I should like to thank Professor Szczepan Szczeniowski for helpful discussions 
concerning the subject of this paper. 
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A METHOD FOR MEASURING THE MOBILITY OF 
IONS IN DIELECTRIC LIQUIDS 
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Physics Department, Academy of Medicine, Gdansk 


(Received July 31, 1958) 


A direct method of measuring the mobility of current carriers in dielectric liquids 
has been developed. The method is based on the measurement of the time of arrival at 
the electrodes of a layer of ions produced by X-rays in a flat ionization chamber 
filled with the liquid under investigation. The velocity of the ions is determined from 
the measurement of the difference in the times of arrival of the ions at the respective 
electrodes for two different initial positions of the ionized layer. The measured values of 
the mobility of the positive and negative ions in hexane are 4.1 10-* cm?/V sec and 
1.3 x 10-* em?/V sec, respectively. 


1. Introduction 


There has been an increased interest in the investigation of the electric properties 
of liquid and solid dielectrics, in particular as regards the dependence of the mechan- 
ism of electric conductivity and breakdown on the chemical structure of the mole- 
cules of the dielectric material (Adamczewski 1958; Lewis 1956, 1957, 1958; Crow et 
al. 1954; Braag et al. 1954; House 1957). These investigations have given rise to_a 
number of problems connected with the determination of basic physical quantities 
(coefficients of mobility, recombination and ion diffusion) describing the mechanism 
of conductivity. i sin : 

The experimental work carried out till now on determining these quantities comes 
chiefly from the years 1932—1937 and must be brought up to date from the view- 
point of recent technical achievements and new problems arising from the recent 
experimental and theoretical work. 

In this connection, the authors, in their investigation on the electric conductivity 
of liquid dielectrics have developed a method for measuring the mobility of ions in 
these liquids. The method described below is a modification of the direct method of 
measurement of current carriers in semiconductors. Its principle is as follows: 

Between the electrodes of an ionization chamber filled with the liquid under 
investigation (Fig. 1), a plane layer of ions is produced by means of X-rays. Next, 
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Fig. 1. Diagram of ionization chamber. A — collecting electrode, B — protective electrode, E — ion 
layer, D — high voltage electrode. 


a constant voltage is applied to the chamber and the time it takes the current carriers 
to reach the respective electrodes is measured. If one repeats the measurement for 


another position of the ionic layer between the electrodes, the mobility is obtained 
from the relation 


As U 
= — LL — l 
EY T E Q) 
where v — is the velocity of a given type of ion; 


u — is the mobility; 

As — is the displacement of the ionized layer; 

At — is.the difference in the times that the ions reach the electrodes in the 
first and second measurement; 

U — is the voltage applied to the chamber; 

d — is the distance between the electrodes. 


A block diagram of the measuring equipment is shown in Fig. 2. 


The closing of key P, switches on the X-ray tube and earths the ionization cham- 
ber K. After a time of about 1 second necessary for producing the layer of ions of sui- 
table density, key P} is opened, which causes the application of the voltage to the 
chamber, triggering of the oscillograph time base and the flip-flop circuit. The use of 
relay P, and a phase splitting circuit at the input makes it possible for the electronic 
circuit to operate under either polarity of the chamber voltage. A delay circuit consis- 
ting of flip-flop circuit I and a differentiating circuit ensures that the voltage on. the 
chamber attains a steady value before the measuring circuit is opened. 

A positive voltage pulse obtained from the delay circuit actuates flip-flop circuit 
II which controls the stage with relay P;. The latter removes the earth from the grid 


‘of the electrometric valve during the time of measurement, this time being adjustable 
within wide limits. 
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~220V =20V 


Fig. 2. Diagram of measuring equipment. S — collimator slit of X-ray beam, K — measuring ionization 
€ 
chamber. 


The preamplifier is a Barth bridge circuit (1934) built with a 954 valve and working 
with a free grid. The control grid current is about 5 x 10-15 A. The output of the 
amplifier is connected to a type OSW-2 oscillograph. The voltage sensitivity of the 
entire measuring circuit is 25 mV/cm. 1 


ton.chamber ionization 


ning measu- vottage on 


voltage ope- 


of time base ring circuit 


triggering 


G 


Fig. 3. Sequence of operations of measuring equipment. 


The sequence of operation is illustrated on the diagram in Fig. 3. The motion of 
the ion layer produced in the liquid by the X-rays in the time t,;+t, begins at the 
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time t, corresponding to the application of the voltage to the ionization chamber. 
As the ion layer approaches the electrodes the potential of the collecting electrode 
increases and attains a constant value at the moment when all the ions arrive at the 
electrodes. The X-ray tube voltage was 60 kV and the current 10 mA. The beam was 
collimated by a slit S of 0.5 mm width and 45 mm in length. 

The ionization chamber (see Fig. 1) was made from brass and the electrodes were 
nickel coated. The distance between electrodes could be changed within the limits of 
0 —15 mm. The surface area of the collecting electrode was 7 cm?. 

Hexane with a conductivity of about 2 x 10-19 Q-1cm! at a field intensity of 
10? V/cm was used for the measurement. 


2. Shape of the voltage pulse 


Let us assume that in a plane ionization chamber (see Fig. 1) at a distance b from 
the collecting electrode a layer of liquid of thickness a, parallel to the electrodes is 
ionized. The charge density c in the layer is constant. 


2 
U,7 4.1410 77 - 


U =1.3x10° -»- 


01 Oe OSE 0.4 E SO. 07 C8. 099 109 11 12 "139 44 


t [sek] 


Fig. 4. Measuring condenser voltage vs time, for motion of negative ions towards the collecting electrode 


The moment, the voltage is applied to the chamber, there arises a current produced 
by the motion of the charge. The current sequence resulting from the motion of one type 
of current carriers of velocity v will consist of two parts. Up to the time the first ions 
of the layer reach the respective electrodes the current i, is constant and has the value 


oav " 


idi ia | BR 
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— 


Similarly, the voltage U, on the condenser C increases according to the relation 


oavt 


As the ions arrive at the electrodes the current decreases linearly and its value at 


LI a . 
the time t = 7 (o al is 


U., = 


ig = dx = I (a — vr) (4) 


and the subsequent increase in voltage is quadratic: 


ov ur? 
U, = Ca c = 5 - (5) 
In order to obtain the resulting waveform for the voltage on the condenser we cal- 


culate with the help of Eqs. (3) and (5) the voltages for each type of ion (taking into 
account their density c, velocity v and the direction of motion) we add up the results. 
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Fig. 5. Measuring condenser voltage vs time, for motion of positive ions towards the collecting electrode. 


Figs. 4 and 5 show the voltage waveform U, for both directions of ionization 
chamber voltage. The values of the parameters used and the polarity of the voltage 
are shown in the figures. It is seen that the displacement of the gap and change of 
polarity of the voltage make possible the measurement of the velocity of each type 


of current carriers separately. 
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3. Results of the measurements and estimation of the accuracy 


Figs. 6 and 7 show the oscillograms obtained under the conditions given for the 
curves of Figs. 4 and 5. In both figures the lower oscillograms correspond to an initial 
position of the ionized layer (4s = 0) and the upper — with the layer displaced. The 
units of the time scale (distance betweeen two dots) in Fig. 6 is 0.02 sec, and in Fig. 7 
— 0.08 sec. The measurements were made at a temperature of 19°C. 


Fig. 6. 


The mean values of the mobility calculated from the oscillograms are 
H+ = 4.1 x 10-4 em?/V sec, 
H- = 1.3 x 107 cm2/V sec, 


The mean square error for ten measurements does not exceed S9 


The-values of mobility for ions in hexane given up till now in the literature, 


12:8 x 10 cm?/V sec and u_ = 4.4 x 10-4 cm?/V sec based on the measurem- 
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ents of Adamczewski! (1937) are smaller than those obtained by the authors. A partic- 


ularly distinct: difference occurs for negative ions. It is probable that this difference 
results from the presence of impurities in the hexane, which are difficult to avoid. 
Reiss (1937) drew attention to the fact that his measurements indicate that the im- 
purities have a particularly large influence on the mobility of the negative ions. This 
problem’ will be the subject of further investigation by the authors. 


Fig. 7. Oscillograms obtained during the motion of the positive ions towards the collecting-electrode 
under the conditions given in Fig. 5. 


As may be seen from Eq. (3) the final increase in waveform voltage is quadratic 
in the time — and this makes it difficult to determine the time ¢, or t5. Since an accur- 
v 


ate measurement of the quantities As, d and U occurring in Eq. (1) does not present 
any difficulty, then the accuracy of the measurement of mobility under the assump- 
y, | 


1 In this work it was also found that positive ions have a mobility 8.2 x 10-4 cm?/V sec. 
= : t . : . H 
However, the number of these ions rapidly was falling off with an ies dn the intensity of 
the M beam. The authors have not observed the presence of such ions. 
E y 
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tion that the ionized layer is parallel to the electrodes, will depend, for a given distance 
between electrodes, on the ratio a/As. The factors governing the thickness of the 
ionized layer are the following: the geometry of the illumination by X-rays, "neutral" 
diffusion, and the effect of the space charge. The last factors were discussed by Reiss 
(1937). It turns out that for the ionization densities below 1071! C/em* used by the 
authors and for ion collecting times of the order of several seconds, the thickness of 
the layer depends on the quality of the collimation of the beam. The (mean) thickness 
is about 1 mm. It should be possible to reduce it to 0.5 mm. 

Thus far in measurements of the mobility of ions in liquids the method of Lan- 
gevin (1903) or a variations of it (Bijl 1912, 1913; Jaffe 1910, 1913; Adamezewski 1934; 


Sato et al. 1956) chiefly have been used. Reiss (1937) developed a method in which the _ 


investigated liquid is moved mechanically. 

In the method of Langevin and in its modification the accuracy is affected, by the 
time measurement, by the stability of the voltage source which produces the ions, and 
the constancy of the sensitivity of the equipment during the runs giving points on one 
curve. Also, the recombination effect, particularly for longer ion collecting times has 
a negative effect on the accuracy of the results. 

In the method presented by the authors, the stability of the ionization source is 
of no significance, and the recombination effect takes place over a relatively short 
time, i. e. only one layer of positive charge coincides with the negative one. In addi- 
tion one can profit from the advantages afforded by oscillographic measurements. 

We should like to express our sincere gratitude to Professor Adamczewski for 
suggesting the subject of this paper and for his frequent valuable advice and aid during 
the measurements. 
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PRELIMINARY -X-RAY INVESTIGATION OF THE STRUCTURE 
THERMISTER MIXTURE OF NICKEL AND MANGANESE OXIDES 
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(Received August 16, 1958) 


The structure of a thermister mixture of NiO and Mn,O, oxides with the percen- 
tage of Mn,O,, changing in the range from 17% to 86.9% was investigated. It was found 
that in the case of all samples there occurs the cubic system of the space group 07 with 
elementary cell dimensions close to those of the spinel of NiMn,O, and changing with 
a change in the percent mixture. In addition to the cubic system of the space group 07 
of the materials investigated there occur other phases different for different percentages 


of Mn,0,. 


The subject of this work was the investigation of the structure of mixtures of 
nickel and manganese oxides according to the relative composition of the mixture. 
These mixtures have a large negative temperature coefficient of resistivity (Schmidt, 
Mikke, Soltys 1955). The samples were prepared by Mr. Schmidt of the Institute of 
Basic Problems of Technique. They were obtained from NiO and Mn,O,. The oxides 
were carefully ground and mixed wet in a ball mill with an addition of caolin and under- 
went heat treatment in a tube furnace at a temperature of 1300°C. 

Photographs of the Debey-Scherrer interference lines of the samples investigated 
were made by means of a chamber of 57.3 mm diameter. Filtered radiation from a tube 
with an iron anode was used. Samples with a thickness from 0.2 to 0.3 mm were inve- 
stigated. In the case of samples which have two phases with similar dimensions of 
their elementary cells, the investigations were repeated using a chamber with a dia- 
meter of 114.6 mm and radiation from a chromium anode in order to separate the coin- 
ciding rings coming from both phases. In the case of the chamber with a diameter of 
114.6 mm, samples with thicknesses from 0.1 to 0.2 mm were used. In the calcula- 
tion of the lattice constant the errors resulting from the shrinkage of the plate and the 
inaccurate determination of the chamber radius were decreased by the Straumanis 
method. A correction for the absorption of the sample was made. The errors caused 
by improper centring of the axis of rotation of the sample in the elementary cell were 
eliminated by the method of least squares. The elementary cell dimensions were de- 
termined to an accuracy of 0.01 À for a cubic system and to an accuracy of 0.02 À 


for a tetragonal system. | 


(199) 


200 Teresa Bedynska 


Fourteen samples of different mixtures were investigated. The crystallographic 
systems obtained and the dimensions of the elementary cells of the individual phases 
for different amounts of Mn,O, by weight are given in the table I. 


Table I 


cubic space group O7 tetragonal 
spinel type 


A plot of the lattice constants occurring for the phase vs the percent mixture 
is given in Fig. 1. 

On all films obtained for different percentages of Mn,O, there occur Debey- 
Scherrer lines of a cubic system with a lattice constant varying from 8.38 A to 8.44 A. 

On the basis of the observed extinctions and reflexes it follows that there exist 
space groups O7, Fd3m. 

For 17% Mn,Og, an accurate measurement of the lattice constant of the cubic 
system of the space group O7 was not made because of the very low line intensity of 
this system. In this case the NiO lines appear very distinctly. 

For 38% Mn,O, in addition to the lines of the cubic system of NiO there appear 
lines of the cubic system of the space group O7 with the lattice constant of the ele- 
mentary cell equal to double the lattice constant of NiO. In the range from 50% to 
86.9% Mn,O; there may be distinguished five regions in which there occur different 
phases. LM 


In the region of 50% (Fig. 2a), 55%, 60%, 65% Mn,O, there occur the Debey- 
Scherrer lines of NiO and lines of the cubic system of the space group O7, with the 


lattice constant of the elementary cell constant within the limits of error and equal to 
8.41 A. é 


At 70% Mn,0, (Fig. 2b) the lattice constant of the cubic system decreases to 
8.38 A and there appear lines of the tetragonal system. The NiO lines do not appear. 
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At 75% and 80% Mn;0, (Fig. 2c) there appear lines of the cubic system of the 
space group O7, lines of the tetragonal system and NiO lines. The lattice constant of 
the cubic system in this region is greater than at 70% Mn$0,. The line intensity of the 
tetragonal system in this region is smaller than at 70% Mn,O, and considerably smaller 
than the intensity of the lines'of the cubic system. 
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Fig. 1 


For 84%, (Fig. 2d), 85.7%, 86% Mn,Oz there occur only lines of the cubic system. 
At 86.6% and 86.9% (Fig. 2e) Mn,O, there occur lines of the cubic system and lines 
of the tetragonal system. m A 

Photographs of the Debey-Scherrer lines for the individual regions ate given 
in Fig. 2. 

In Fig. 3 are shown the regions of different phases on the background of the 
curve of resistance vs percent mixture. The plot of resistance vs the percent mixture 
was taken from the work: done by Schmidt, Mikke, Soltys (1955). 
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According to data from the literature, Mn,O, forms from Mn,O, at a temperature 
of 1080°C (Kage, Lavy 1956). Mn4O, has two crystallographic forms: tetragonal 
(Wyskoff 1951) and cubic (Murdie, Sullivan, Maurer 1950). The cubic form has 
a structure of the spinel type with the length of the edge of the elementary cell equal 


— ——». 


Fig.2a,b,c,d,e 


p 


to 8.64 À (Murdie, Sullivan, Maurer 1950). The tetragonal form has a deformed spinel 
structure. It can be described by means of a tetragonal system with elementary cell 
dimensions a = 8.14 A, c= 9.42 A (Wyskoff 1951). The tetragonal variant forms - 
at a temperature of 1080°C, the cubic one forms from the tetragonal form at a ms | 
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perature of 1170°C. A NiMn,O, spinel with an elementary cell edge of 8.38 A 
(Frevel 1942) forms from a mixture of nickel and manganese oxides. 
The cubic system appearing on all plates has elementary cell dimensions close 


to the dimensions of the elementary cell of NiMn,0,. 

The amounts reqiured to produce NiMn,O, are 67.89% Mn,O, and 32.18% NiO, 
which corresponds to the molar ratio 1:1. 

At 70% Mn,0, the lattice constant of the cubic phase is equal to the lattice 
constant of the spinel of NiMn,O 
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Fig. 3 


In the region 50%, 55%, 60%, 65% Mn,O, there appears a cuba system of ue 
spinel type with a lattice constant a little larger than for pure NiMn,O,. This is Br e 
the spinel of NiMn,O, with the lattice constant changed de cue i ; 
It is possible that the cubic system observed in the region 75%, 80%, 84 vos i Ve: 
86%, 86.3%, 86.9% Mn,O3 with a lattice constant greater than the me nd 
of pure spinel of NiMn;O, is a solid solution of NiMn,O, En the wen BEE: 
The tetragonal phase appearing at 70%, 75%, 80% 86.3%, 86.9% x A $ = SS 
owing to the deformation of the tetragonal form of Mn30, by the a s 2 à ). 

Conclusions: 1) For all samples investigated there appear the a e 1e cubic 
system of the space group Oj with elementary cell dimensions close to the dimensions 


of the elementary cell of spinel of NiMn,O,. 
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2) The lattice constant of the cubic system changes with a change in the percent 
mixture. | ree: 
3) In addition to the eubic system of the space group O% there appear in the 
materials investigated different phases dependent on the percent mixture. 

4) For samples showing a minimum resistance there occur two phases: cubic 
and tetragonal. The lattice constant of the cubic system attains a minimum. 

I wish to thank Professor L. Sosnowski, Mr. J. Auleytner for discussions of the 
results and constant aid, and Miss. Z. Wójcik for help in the experimental work and 
in the calculations. 
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This work deals with an investigation of surfaces of semi-conductor monocrystals 
and samples of polycrystalline indium antimonide by means of etching with a suitable 
etching mixture. The etched figures obtained on the surfaces of different crystallographic 
orientation had the shape of furrows changing in width and shape. The patterns obtained 
on the surfaces are not the result of the existence of a superstructure. Their occurrence 
can be explained possibly by a tendency of the boundary of the phases of the crystal 
surface and the etching medium to attain minimum energy. 


The work described below is the result of an investigation of cut and ground 
surfaces of the semi-conductor indium antimonide by means of chemical etching with 
a chemical agent affording fine etching. Mono- and polycrystalline samples were 
used. 

The action of suitable etching agents on metals and alloys can cause the formation 
on their surfaces of patterns in the form of furrows which change in width and shape. 
Surfaces of this type can be obtained both by electrolytic and by chemical etching. 

K. F. Hulme (1955) electro-etched zinc monocrystals and obtained on the surfaces 
of the samples patterns in the form of furrows. Similar figures were observed also on 
crystals of an aluminium-manganese system by T. Imura and T. Suzuki (1955). 

In this work indium antimonide monocrystals of a purity corresponding to 10!® 
holes/cm? (type p) were studied. / i 

The samples were prepared in the following manner. A suitably oriented cut 
surface was ground with sand paper (with carborundum grains) of numbers 220, 320, 
400, 600 successively, then polished with felt and etched by an atching mixture of the 
following composition: 

j 1) HF, 5 parts 
2) HNO, 5 parts 
3) CH,COOH 3 paris (saturated with quinine) 
4) CHOH 2 parts 
: 5) H,O, 6 parts. 
(205) 
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The etching time was about 6 hours at a temperature of 15°C. The etched sur- 
faces were rinsed in methyl alcohol and dried. 

The orientation of the various surfaces was established by means of reverse 
lauegrams. A point source of X-rays with a focus diameter of about 30 u (Auleytner 
1957) was used. : 
` The prepared surfaces were observed under a metallographic microscope with 
a magnification of x 250 to x 1200. The images obtained showed differently oriented 


Fig. 1 


surfaces with furrows whose shape and width on the same surface changed either 
continuously or in jumps (furrow density from 10 to 40 per 100 u (Fig. 12). 


In addition to the above there were observed local increases in furrow density 
in some regions and decreases in density in others (Fig. 3). 


The surfaces of polycrystalline sample cross sections were covered with furrows 
of different orientation and density (Fig. 4). 


Also observed were inclusion strips on which the furrows were. arranged per- 
pendicular to their boundaries and marked by a good deal of regularity (Fig. 5). 
It was noted that on wider strips the degree of ordering in the furrow system was 
smaller than on the narrow ones. As shown by the reverse photographs of the Laue 
method, the strip inclusions were twins. 


In some regions the furrow distribution was analogous to the equipotential lines 
about two electrostatic charges of the same sign, where the distances between centres 
of the disturbances were about 250 u (Fig. 6). The density of the furrows was measured 
along the normal to the line joining the two disturbances. A plot of the furrow density 


vs the distance from the centre of the line j joining both VENDRE was of a hyper- 
bolic character Os. a 
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Fig. 2a 


Fig. 2b 


It was noted that in the polycrystalline substance the furrows on some of the 
edges of the grain penetrate into the neighbouring crystalite to a depth of several 
microns, while others did not (Fig. 8). In some regions there appeared systems of 
furrows making an angle of nearly 60° with one another (Fig. 9). The X-ray photo- 


graphs did not show a difference in the orientation of the regions on both sides of 


the boundary formed by the two systems of furrows. 
On surfaces close to the crystallographic plane (111) (at an angle of about 10°) 
| there appeared pits several microns deep and covered with concentric furrows. Their 
symmetry corresponded to the symmetry of the crystallographic direction of the given 


plane (Fig. 10). > 


B 
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Fig. 3 


Fig. 4 


It was found that by using an etching agent of another composition the shape 
of the etched figures changed (Auleytner, Kolakowski 1958). 

The reason for the occurrence of surface figures of this type, as given by various 
authors, are different for different materials. While T. Imura and T. Tomizuki believe 
that for the aluminium-manganese system they may reflect an existing superstructure, 
K. F. Hulme reported that for electro-etched zinc the density and shape of the furrows 
changes with the intensity of the current used for the electro-etching. 

Since in the case of the indium antimonide crystals under investigation here no 
superstructure was formed and no current was used during etching, it seems that the 
causes of the formation of figures of this type must be explained in a different way. 
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Fig. Sa 


Fig. 5b 


+ 


The formation of furrows on the etched surface can be explained by a tendency 
of the boundary of the phases — crystal surface — etching agent — to attain minimum 
energy. The etching process takes place in the case of indium antimonide in two stages: 


1) formation of a layer of oxides, 

2) dissolving of the oxides that have been formed. On the different crystallo- 
graphic walls that form in the first stages of the etching the oxide films have different 
properties because they from on different bases. In the case of an etching agent of 
appropriate composition some molecules can be adsorbed on some specific surfaces 


and to a smaller or larger extent impede the etching process. 
p ' 
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OI E 
In the given case it is most probable that the quinine molecules, being one of 

the components of the etching solution, are adsorbed on the surfaces having a specific 

amount of energy, at lower values of energy the quinine molecules are weakly adsorbed, 

in which case they do not separate the crystal surface from the etching agent and do 


Fig. 6 


not impede the etching process. The occurrence of a surface of too great an energy 
is not possible, since the adsorbed molecules do not decrease the surface energy 
sufficiently enough to impede the etching and the region close to the surface tends 
to the most stable state corresponding to a minimum surface energy (Fig. 11). 


` 
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Hence the impeding of the etching as a result of the adsorption will occur on 
surfaces of constant energy and the furrows denote the isoenergetic (surfaces) in the 
crystal. 


the furrows are not necessarily the crystallographic planes. The furrow density depends 


If the crystal does not have any stresses the furrows are restricted to the crystallo- - 
graphic planes. In case there exist internal stresses in the crystal the surfaces limiting. 


PU We 
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Fig. 8 


on the amount of etching nuclei which form during the etching process and on the 
ratio of the velocities of their lengthwise to crosswise growth. If this ratio is small 
the furrows are very wide (15, Fig. 12) and if it is large the furrows are narrow (2, 
Figs. 1, 5). 

_ With an increase in the density of the nuclei the shape of the furrows becomes 
more regular. A factor modulating the shape of the furrows are the different stresses 
on different regions of the crystal surface. They arise during the growth of the crystal 
through: the joining of atoms and also of one- and two-dimensional crystals, and three- 


dimensional nuclei. gj 
s 
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As a result of their imperfect attachment to the growing crystal there arise stresses 
which reach far beyond the deformed region. They appear during the etching in the 
form of more dense and less dense furrows. This is illustrated schematically by Figs. 
13a, 13b, 13c, and 12. 

These stresses can have a dislocation character, especially on the boundaries 
where the two-dimensional crystals adhere. 

The observed stress regions extend up to about 200 u (observations of Bond and 
Andrus 1955) for edge dislocations of about 100 u. 


II. Quinine molecules are strongly adsorbed and, as a result, the energy of the boundary between 
the phases decreases — the etching being considerably impeded. 

m. Quinine molecules are’ strongly adsorbed, but the energy of the boundary of the phases — 
wall 3 — etching solution — is so great that the adsorbed quinine molecules cannot be completely satu- 
rated and consequently the etching proceeds until the wall disappears. ae 


i 
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Fig. 12 


a b c 
Fig. 13 abc > 


In certain regions there appeared furrows forming a closed polygon (Fig. 14). 
Their shape is probably a reflection of the stresses arising in the crystal lattice 


around the dislocation. 


Conclusions 
1) The etched patterns on the surfaces of indium antimonide monocrystals in 
the form of furrows are not closely connected with the crystallographic planes and 
they appear on surfaces of different orientation; they are probably a reflection of the 
stresses arising in the crystal during its growth and denote the. surfaces (directions) 


of the same energy properties. 
2) The shape of the furrows and the space between them change in different 


< 


y 
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regions of the same crystal surface. The distance between neighbouring furrows varies 
from 2 u to 15 u. The regions of disturbances in the parallelness of the furrows extended 
from about 15 to 200 microns. 

3) Some regions were covered with furrows whose density distribution resemble 
the density of equipotential lines about two electrostatic charges of the same sign, 
but the change in density along the normal through the centre of the line joining 
the centres of the disturbance was of a hyperbolic character. 


Fig. 14 


4) It was found that on twin inclusions in the shape of narrow strips the furrows 
are very regularly situated and that with a decrease in the width of the twin the regu- 
` larity and density of the furrows increase. 

The author would like to thank Professor L. Sosnowski for discussion of the work 
and for his valuable comments, Mr. J. Auleytner for his advice and remarks, Mr. W. Gi- 
riat for supplying the indium antimonide monocrystals grown by him, and Mr. J. Rautu- 
szkiewicz for determining the purity of the crystals. 
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In order to investigate the spin waves in ferrimagnetic translation lattices, a new 
classification of superstructures is given, introducing the definitions of natural order 
and regular orders. By means of a semiclassical approach the general dispersion formula 
for spin waves frequency is derived, taking into account interactions with all neighbours 
in any regular spin order of translation lattice with two kind of spins and comprising also 
ferromagnetism and antiferromagnetism. In the special cases the derived formula agrees 
entirely with those of other papers, based on semiclassical thcory. The kinematical method 
of Keffer, Kaplan and Yafet is generalized and his equivalence with formal calculations 


is shown. 
The picture of spin precession waves in ferrimagnetic translation lattices with two 


kinds of spins is investigated. 

The analysis of non-regular superstructures leads to an interesting consequence 
for antiferromagnetics. It appears that simple antiferromagnetic sublattices in body 
centered non-regular superstructures must be mutually uncoupled. This result is confirmed 


for the crystal of MnO,. ; 4 


I. Introduction 
The purpose of the present paper is a discussion of spin waves in various types of 
ferrimagnetic superstructure in translation lattices with two kinds of spins and a deri- 
vation of the general formula for the spin waves dispersion, which should be valid 
for a large class of ferrimagnetic orders, called here “regular orders”. We use for this 
purpose the semiclassical approach, since it is a very convenient method for investi- 
gation of spin waves propagation. Moreover, this approach leads to the same expres- 
sions for frequency spectrum as the quantum mechanical theory of spin waves in the 
Holstein-Primakoff approximation, as pointed out by Kubo (Nagamiya et al. 1955, 

page 102) for the antiferromagnetic case. | 

(215) 
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1. The semiclassical approach 


The semiclassical approach (Heller, Kramers 1934) consists in treating of spins 
as classical vectors, coupled by the exchange energy 


H IS 2] Ss (1) 


m,n 
the Jinn being the exchange integral for interaction between spin vectors AS,, and AS, 
h being Planck's constant divided by 2%. Such an approach leads to the precession 
wave picture of spin wave, because from (1) we obtain the following equation of mo- 
tion of the m-th spin (Dàring 1948, Herring, Kittel 1951): 


; 2 
SESS Sa DS 2) 


using the abbreviation: J, = J mm +r The summation extends over all lattice vectors 
r, stretching from the site m and being included within a considered ferromagnetic 
domain. The symbol P = yy, (Ho + HL,) contains the external magnetic field Hy as 
well as the effective field H , of magnetic anisotropy, y standing for magnetomechanical 
ratio and u, being the permeability of vacuum (we are using Giorgi's units system). 


It is easily to verify, that the vector 


Sm = R cos (km + ot) j, + R sin (km + ot) j, + 5* j, (3) 
is a solution of (2), when the relation occurs 
w = —P --25h-1S* > iR (cos kr — 1) . 4 (4) 


provided, that the lattice is translational. Here the j,, Jy j, are the versors of carthesian 
right system of-coordinates axes with axis z parallel to P. One may see from (3), that all 
vecto:s S,, precess about P-axis with angular velocity œw and with radius R of preces- 
sion. Two spins whose positions differ by * have a phase difference O = kr in their 
precession. 

In a similar way one may show, that in ferrimagnetics the motion of spins 
forms also a precession wave with fixed phase-relationship between different sublattices. 
Using the solution of form (3) one must consider that R and S* both depend on lattice 
site 3? ; moreover, one must introduce also a phase difference a: The dependence of R 
upon m appears even for antiferromagnetics, as demonstrated by Keffer et al. (1953). 
Substitution of such a modified expression into (2) shows, that under some conditions 
for w (equivalent with dispersion law), this is a correct solution of eq. (2) in the case 
of ferrimagnetism. Thus, for ferrimagnetics we get the same picture of spin waves, 
as for ferromagnetics. We shall proceed to a more detailed discussion of this spin 
wave picture in subsequent paragraphs. 


2. Two ways of calculation 


There are two ways of treating of spin waves, based upon the semiclassical equa- 
tion of spin motion. The first way treats this equation formally, as shown for ferro- 
magnetics in (3) and XR Substituting the guessed (e. g. Dóring 1948. or Cofta 1956) 
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or calculated (Nagamiya et al. 1955, page 98) wave precessional expression into the 
equation of motion, one obtains certain conditions (for the direction of precession 
axis, for the allowed phase differences and for the relations between frequency and propa- 
gation vector, i. e. for dispersion law), under which this equation is satisfied. This way 
was used by many authors for some special cases of ferrimagnetism, namely by Kaplan 
(1952) for ferrites, by Kouvel and Brooks (1954) for simple cubic lattice, by Kouvel 
(1955) for magnetite and by Kubo (Nagamiya et al. 1955) for natural orders. All this 
authors write the equations of motion separately for each sublattice (likewise Cofta 
(1958 a) for the case of lattices with three kinds of spins), but that is not necessary, 
since the equation (2) comprises all cases. 

The second way is the well-known kinematical approach of Keffer, Kaplan and 
Yafet (Keffer et al. 1953), generalized in the present paper for any translation lattice. 
It uses the precession equation of a single spin: 


S—-oxS (5) 
which, according to the character of solutions (3), must be fulfilled by each spin of the 
lattice. Let us introduce a moving carthesian coordinates system revolving together 


with spin, directed in the manner shown at 
fig. 1. Then, we may write: 


= wt R (6) 
It is easy to calculate from (6) the w* for a given 


sublattice, since one can find S" from the 
appropriate equation of motion. If the same 
angular velocity for both sublattices is required, 
we get the conditions for œw, leading to disper- 
sion law. 

Hence, the kinematical approach allows us 


to derive dispersion laws for spin waves. Kef- 
fer, Kaplan and Yafet (l. c.) have applied this Fig. 1. ~ 
approach to explain the mechanism of spin wave 

propagation in antiferromagnetic linear chain. It is used also by Cofta (1956) for 
ferrimagnetic linear chain. : 


3. Classification of spin orders 

In ferrimagnetic lattices with two kinds of spins the most simple order is the one 
in which any spin of first kind is surrounded by spins of second kind = and vice 
versa. Let us call this.order a “natural order”. 

We shall distinguish now a larger class of orders. Let us call a “regular order” 
each superstructure, in which each lattice point is a center of symmetry (infinite lattice 
is assumed). That means, that the two neighbours in sites r and —7 relative to an arbi- 
trary spin belong both to the same sublattice. Note, that the natural order in ferrima- 


_ gnetic translation lattice is evidently regular. 
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Let us introduce a condition 
er er TT (7) 
If this condition is assumed, the abbreviation J,,,,,, = J, used in (2) for ferro- 
magnetic case is correct also for any regular ferrimagnetic order. 


II. The dispersion law for regular spin orders 


We shall derive the general dispersion law valid for any translation lattice with 
two kind of spins, ordered in the regular manner!. We base on the equation of motion 
of the m-th spin 


Sm = Sin x Uu oe 2h > Jy "— (8) 


with a vector P. dependent on lattice site. Moreover, we impose the Born-Karmán's 
m p P 
boundary condition, called also “periodicity condition". 


1. The formal derivation 


According to paragraph 1.1 we use the precession-wave solution, writing it rather 
in the convenient complex form 


Scu ug) c js (9) 
Substituting (9) into (8) and comparing appropriate vector components, we obtain 


h (o + PZ) Ry, em = — 225: y p J, Rn em + SE 2X], R, ,, 69" ** e". (10) 


m+r 


Now, we put for the first sublattice of a regular superstructure 


Rea Ra, = 005, sm 
z z z z | (11a) 
PE P= ys (Ho + Hy) 
and for the second one 
Rm = Rs, An = a, Im = B, | (11b) 
Pm = Q* = yuo (Hg + H3) 
the H, being the anisotropy field for i-th sublattice. Using abbreviations: 
+ + 
Ji (k) = 251 X J (er — 1) = 2n > J, (cos kr — 1) 
r T 
fa (k) = 2h71 X} J, er = 253 X J, coskr (12) 


x-— —2h1 Y jJ, : 


where i symbolizes a summation over all lattice points belonging to the same sublat- 


1 That is not true when not all points of the same sublattice are geometrically equivalent, 
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tice as the considered point m and D symbolizes a summation over the other lattice 
pointe, we obtain the simultaneous equations 


R, (A*f, + B*x — P* —9) + Rye A*f, = O 
13 
R, Bf, + Raes (BA + 4*y — Q* —w) =O 23 
The secular equation of (13) leads to the two values of angular velocity 
1 z z z z 
Ed BU M DE 
; (14) 


+ 3 (AF-B) (ef) + (P'O) + 4A BD 


corresponding to the two branches of frequency spectrum. Note, that the secular 
equation contains not the phase factor exp (ia); this means that the dispersion law (14) 
does not depend upon phase difference a. 


2. The kinematical derivation 


In the case of regular superstructures, fulfilling the condition (7), we can write 
also the eq. (8) in the form 


h S, = S Xb Py +2 Ip Wanye (15) 


Each vector 
= S% x Sie, aeg ») 


is always perpendicular to the plane, formed by S, and œ; hence, in the system shown 
at fig. l, we have 


W»,-—2(StR 


mio COS Rh Sta RL) (16) 
. puting here both a,, = 0, since dispersion law does not depend upon phase difference. 
We take also into account, that as can be easily shown œ must be parallel to P. Then, 


assuming the same angular velocity for both sublattices 


wi = 03 = (17) 
we obtain from (6), (15), (16) and (17) the expression : 
ho Rim 2a) Stt epcos kr. ei Hop (18) 


from which the dispersion law can be derived in a way similar to the formal one as 
shown above. 

Thus, the generalized kinematical approach of Keffer, Kaplan and Yafet was 
applied here for derivation of dispersion law for arbitrary ferrimagnetic translation 
lattice, ordered in a regular manner. The approximations for small k or small R were 
not necessary. 

The expression (18) is identical with (10) for a,, = 0, since the sums arc the same 
for regular orders (see (12)). Then it leads to the same result (14). We see, that both 
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the ways, the formal and the kinematical one, are equivalent as methods of derivation 
of dispersion law. It may be easily shown that the (10) and (18) are also identical for 
any phase factor exp (ia). 


3. Special cases 


We shall consider now some special cases. of dispersion law (14) in absence of 
an external field (#, = 0). For this purpose, we introduce the symbols: 


|47| E S |B?| = Se, || = H; 
(Note, that all these symbols represent the z-components.) 
a) Parallel spins 
In this case we have 
A= Su. Bi 53, (Ph — vito O — YHH. 
Thus, the formula (14) takes the form 


w= 5 (8, +5) % +A) Hr (Hy + H9) + 
£5 (US, — 5) & A) — yo He HS 9) 


When in a natural order we restrict ourselves to the nearest neighbours interaction 
only (f, = 0), then, neglecting the anisotropy field, we find 


l ~ N 2 2 i 
nr (Si Se) + p Ue (Sp — S1)? + 4 5,5, fiy (20) 


in agreement with expressions (15.5) obtained by Kubo in the paper of Nagamiya et 
al. (1955). In the case of simple cubic lattice with lattice constant a we obtain from (20) 
an approximate formula for small k (we expand the square root): 


EJ Ry oder S15; 
o c 6I (Sy +S) (1-1) AJ ZT ate (21) 


One of two branches (21) has been obtained previously by Cofta (1956) for linear i 


ferrimagnetic chain. 
When 5, = S, = S* and Q7 = P7, we must put consequently f, = O and x = O. 
We have now the case of ferromagnetism for which we obtain from (19) the relation (4). 
b) Antiparallel spins 
In this case we have | 
Pe à B= —S,, P? = —yu, OF = yuH, 
Then, (14) changes into 


o = $5 — S) (A +) + yuo (Hy — Hy} + 


i | 
+ y US +5) 05 — 9) + yuo (H + HOP — 45,5, 2) Q2). 


a 
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Taking into account the interactions between nearest neighbours only and neglecting 
the anisotropy fields, we find from (22) for a natural order the formula 


"rl 
(Sy = 5) Ez {x (Sy + S)? — 45,5, £334 (23) 


which agrees with formula (15.10) derived by Kubo in the paper of Nagamiya et al. 
(1955). On the other hand, for a simple cubic lattice ordered in natural manner, ta- 
king into account the external and anisotropy fields, we obtain the formula! (2.7) 
of Kouvel and Brooks (1954). In the same case, for small k the formula (23) takes 
the approximate form: 


o 62 8-5) £15, 812 2 Poe (24) 
3 |S; — S, 
The antiferromagnetic case is excluded here, because we cannot expand the square 
"root in this case. 
In the antiferromagnetic case we have always S, 4 Sa for k O because R, differs 
from R, (see Keffer et al. 1953). Nevertheless the differences |S,—.S,| are very small 
and we may put 5, = S, = S and H, = H, = H}. Then from (22) we obtain 


w = + VIS (fi —») + yuo Hal? — Sfr (25) 
In the case of natural order, the expression (25) agrees entirely with formula (15.11) 
of Kubo (Nagamiya et al. 1955). For very long waves one can put k = O whence follows 
fy =O and f; = — x. Then (25) leads to the well-known formula (Kittel 1951, Naga- 


miya 1951) for resonance frequency 
lo] = ly| Ho V H4 (H4 + 2Hg) (26) 


where yu, Hg = — x5: 

The formula (25) applied to several antiferromagnetic superstructures considered 
in a Letter to the Editor by Cofta (1957) leads in approximation for small k to all the 
four relations, obtained in the mentioned publication (when one neglects the anise- 
tropy fields and all interactions except these between nearest neighbours). 


4. Some remarks 

For the derivation of dispersion formula (14) we have used one approximation : 
only, i. e. the semiclassical approximation. 

The formula (14) has been obtained also by Cofta (1958 a) as a special case by 
consideration of ferrimagnetic lattices with three kinds of spins (also by means of 
semiclassical: approach). 

The considered formula takes into account interactions between all pairs of locali- 
zed atomic spins. In some cases it is necessary to consider interactions between second 
(i. e. next nearest) neighbours (superexchange). We have such a situation when investi- 


1 We take into account the discrepancy between y-factors in § II. 4 of the present paper. 
¥ 
u 
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gating anisotropic spin wave dispersion in non-natural superstructures (the “super- 
structural anisotropy“, see Cofta 1957 and 1958 b). In such superstructures there 
exists probably always the superexchange between second neighbours. The formula 
(14) has a form very convenient for detailed calculations in the case of dispersion 
anisotropy. Nevertheless, we can write (14) also in another form, using Luttinger’s 
vectors (Ziman 1952). 

For each regular superstructure! there exists the Luttinger’s vector w with the 
following property: the number exp (iwr) is equal to 1, when the vector r connects 
two lattice points, belonging to the same sublattice, and it is equal to —1, when r 
connects.two points belonging to different sublattices. Using Luttinger's vector, we 
can write in (14): 


x f ht» J, A H er) cos kr — 2} 
x — fí— h y J, {2er — (1 + e") cos kr} 
Ufo bh J, (1 — e) cos kr 


We can take also into account in (14) the spin-orbit coupling effects in the manner 
shown by Wangsness (1956), introducing different magnetomechanical ratios y, and y; 
for different sublattices. Using this approach we have to write in (14): 


i P = y, po (H, + H,), Q = yz uo (Ho + M) 
(Compare with the paper of Kouvel and Brooks, 1954). 


III. The picture of spin waves in lattices with regular spin orders 


We shall consider now some kinematical properties of precession waves in trans- 
. lation lattices with a regular superstructure. For simplicity, we may put P = Q = 0. 
Introducing the abbreviations 


Ar sr Den to HE US 
1 ; A* f, = b, Bf, + A* x= d, 
we write the equations (13) in the form: 


(a — o) R, + bR, e* = \ 


cR, + (d— a} N (27) 
whence 
Rast sty 0 

Wie em: (28) 


Since all coeficients a, b, c, d, the angular velocity w and the radii R} and R, are real 
the factor exp (ia) is also real (a = 0 or a= 2). 

Let us denote by c, the value of angular velocity, corresponding to the upper 
sign in (14), and by wœ the other one. 


1 Sce note below the page 218. 
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l. The phase difference 
From (28) we have 


eia Rı see 
R, b 
In the case of parallel spins J > 0; hence 
w— a w — a 
b >O and b <O 


Thus for the first solution (i. e. for the upper sign in (14)) we have a motion without 
phase difference and for the second one with the phase difference a — zr. 
In the case of antiparallel spins it is J < 0 and we have 


w—a 


b «0 


for both solutions. Thus, the spins in different sublattices differ in phase by z. 
2. The inclination of spin vectors 
Let us denote the inclination angles of spin vectors to the z — axis by 9, and ù, 

respectively, as shown in fig. 2. We have then 

R, = Asind,, A* = A cos à, 

R, = Basin 8, B* = + B cos 0,, 
. the upper sign being valid for parallel spins and the lower one for antiparallel. Substi- 
tuting these relations into equations (13) (with P7 = Q* = 0) we find easily, that the 


equality 9, = 2, occurs only in the case of ferro- 
magnetism, or when 


E Z 
Be O (29) 2 
This condition is satisfied e. g. for k — O. It is to be #7 
noted, that (29) is not a sufficient condition. B 
We shall investigate now, which angular in- Y 
clinations correspond to the values w, and w of 
angular velocity. For this purpose we will consider N 
the expression 7 
sind, _ AR, _Alo-a loa (30) 
sind, GBR BAS bi fe a b 
The last approximative equality follows from the Fig. 2 ab 
relation B |b| ~ A |c| valid for small inclinations ®. 
For small values of k we can write approximately (see Appendix A): 
o a (e- b le B) r (31) 


the upper sign being valid for œw, and the lower one for wg. 
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a) Parallel spins 
In this case we have c > 0 and b > 0. The angular velocity w, corresponds to the 
wave with 94/9, = l, since from (31) we obtain 
lo — a | e e| 
In the limit of infinitely long waves (k = 0) we have exactly |w; — a| = |c| and conse- 
quently 9, = d,. 
The angular velocity œ, corresponds to the precession with R, ~ R,, since 
from (31) follows 
Io; — a| ~|b| 
As we see from (30), in this case the shorter vector has always a greater inclination, 


in agreement with the geometrical interpretation (fig. 3). For k = 0 we have exactly 


|o, — a| = |b], that is Ry = Ry. 


Z 
Rs 


>, 


R2 


% 


Fig. 3 


b) Antiparallel spins 


Assuming A* > 0, B* — 0 and moreover |A*| > |B*|, we have (since J < 0) 
b <0, c>0 and |b| > |c]. Consequently b +c < 0 and 


o-a- m {o-b = (c4-B) (32) 


The angular velocity w, corresponds to the precession with nearly equal radii Rover Ro, 
since from (32) follows 

Io; — al =~ |b| 
For k = 0 this equality is strongly satisfied (fig. 4). All details are analogous to the 
case = wy for parallel spins. 

The angular velocity c corresponds to the motion with 9, = Da, as follows from 
(32). For k = 0 we have exactly 0, = 0,, for kÆ 0 however a small difference Da — 9, 
occurs, the sign of which depends on the values of A and B, on the type of superstruc- 
ture and on wave direction. 
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IV. The spin-precession waves in lattices with non-regular spin orders 


It is an interesting question to analyse the semiclassical picture of spin waves 
in such superstructures of ferrimagnetic translation lattices, which are not regular. 
l. General remarks 

The main objective of such a analysis is the non-regular ferrimagnetic linear 
chain of type: 

ACoA eA, BEND 

D BD 8 (9 @ (Q2 (3 @ (33) 
This chain is the simplest example of a non-regular order and moreover it forms the 
basic element of all three-dimensional non-regular superstructures. It consists of 4 
ferromagnetic “sublattices“, marked in (33) by numbers in parentheses. All points 
of the same sublattice are physically equivalent, i. e. have identical neighbourhoods. 
Two sublattices are occupied by A spins and two by B spins. Nevertheless, the preces- 
sion radii are in principle different for each sublattice. 

Such an ordering of spins can form an obstacle for the propagation of considered 
precession waves. We can show this using equation of spin motion for any linear chain 
with nearest neighbours interactions only: 


hS,—2]8 SIS. 348, D (34) 
When the spins S,,_ , and S,, , , belong to different sublattices (which occurs in a non- 
. regular chain), their sum S = S,,_, + S, , , generally does not lie in the plane JI 


formed by vectors S, and œ. Then, S,, is not perpendicular to /7 and the precession 
wave can not exist. This situation is realized in the special case of antiferromagnetism. 
The sum S has always a non vanishing component (R, ., + R,,,):sin ak, 
perpendicular to JJ. In ferrimagnetic case (|S,,_ | |S,, , |) however one can 
choose all radii so, that the sum S lies in the plane IT. 


2. Analysis of linear chain 
Substituting a solution of type (9) into (34) we obtain following equation for 
arbitrary linear chain with constant distance a between neighbours 


hz, Ry etm = Rye (Sp -1 + Sma) + 

ae SERE eltm-1 er ak RE em +1 pu) 

When the point m belongs to the 1-th sublattice, we have to put 
= RR a it Aan, 


m 


(85) 


and obtain thus a system of four simultaneous equations for four unknowns R, exp (aj). 
Their secular equation gives 4 values of angular velocity, namely: 


o,a = J LV + [V° 20 + 2D) 


36 
or J =V - [V* - 2U + 2/DY) sd 
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where 
D(k) = U? — AW(k) 
U = SESE 45385 4 SISE + S151 
V = S + S3 + S3 t Si 
W = S S Se 94: 4 ein? 2ak 
Approximately one can put 
Si =~ S3 = A’ 
$32 94 = B* (37) 
Then 
= (A* + B*)* — 16(A*B’)? sin? 2ak 


This expression is positive for 
(A* + B*)? 


[sin 2ak| < 
4 A*B*| 


(38) 
In the case of parallel spins this equality is satisfied for all k. In the case of antiparallel 
spins it is satisfied for k sufficiently small. In ferrimagnetics with antiparallel spins 
we have always A? + B* > 1, therefore within limits of validity of spin waves theory 
the condition (38) is satisfied. In antiferromagnetics however A* + B*— 1/N (N — num- 
ber of points in chain, see Keffer er al. 1953), hence 


Isin 2ak| < (2S*N)-? < 1/N? 


and |ak| < 1/(2N?). This inequality is valid only for k = 0, because the smallest non 
vanishing value of |ak| is 2zt/N (periodicity condition). 

By the approximation’ (37) we have V?—2U = 2(A4* + B7}. Since 
D < (A* + B?)*, the expression V? — 2U + D* is always positive.’ 

Thus, the results for non-regular ferrimagnetic chain, obtained by semiclassical 
approach, are following. In ferrimagnetic chains (A + B) the spin waves can propa. 
gate without restrictions. In antiferromagnetic chains however, the spin-precession 
waves with |k| > 0 cannot propagate, since their frequency had to be complex. 
. A precession-wave motion in such a chain is possible only inside each of two identical 
antiferromagnetic sublattices separately. (In (33) the one antiferromagnetic sublattice is 
the (1) + (3) sublattice, and the other one is the (2) + (4) sublattice). y such a case 
of uncoupled sublattices one would have to put J, = 0, Dant 0 J, = Jm,m +1) Which 
would be entirely irreal for a linear chain. 


3. Analysis of some three-dimensional lattices 


In every traslation lattice with non-regular superstructure the chains (33) of 
nearest neighbourhood exist. When not all lines of nearest neighbourhood are non- 
-regular, then the restictions similar to these from previous paragraphe concern only 
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some directions. Instead, when all nearest neighbours lines ar non-regular (the "com- 
plete non-regularity”), the situation is identical in principle with that for a non-regular 
chain. An example of such a superstructure of body centered cubic lattice is shown 
in fig. 5. This lattice consists of two antiferromagnetic simple cubic lattices. Each of 
these simple lattices has a natural order. All four nearest neighbours directions (along 
four body diagonals) form non-regular chains. 

The calculations analogous to that made for a non-regular chain lead to the con- 
clusion, that in an antiferromagnetic bcc lattice with "completely non-regular” order 
the precession waves cannot travel through nearest neighbours; nevertheless, they 
can travel through second neighbours (next-nearest neighbours), when J, = 0 but 
Ja 75 0. That means, that the two simple lattices of which the bcc lattice consists 
are not mutually coupled, and consequently the precession axes in different simple 
lattices can have different directions. 

The situation in orthorhombic body centered lattice is identical. 

In the nature there exists one antiferromagnetic lattice with the above considered 
"completely non-regular" superstructure. It is the lattice of Mn ions in crystal MnO,. 


Fig. 6 > 


It is a tetragonal body centered lattice. According to the neutronographic measure- 
ments of Erickson (see Nagamiya et al. 1955 page 7) each of component simple lattices 
is here antiferromagnetic and has a natural superstructure, but the direction of spins 
(identical with their precession axis) in one component lattice is perpendicular to the 
spin direction in the other one (see fig. 6). Thus, both the simple lattices are not 
coupled mutually consistently with our conclusions. 
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Appendix A: The derivation of equality (31) 


In the paragraphe III. 2 we have used the approximative equality (31). In order 
to show the validity of (31) for small values of k, we must take into considerations, 


that 


(except in ferromagnetic case who x — 0), and 
[fa + *| < Ixl 
since all cosines in these expressions are very near to 1. Thus, we find easily, that 
me 7-6. and d= —b 
and consequently 
(d — a)? + 4be z (c + b)? 
Thus, we obtain (31): 


0 —a = > ((d— a) + y (d — a)? + 4bc } => {eb + |e + dj} 


Appendix B: The elimination of precession radii in the manner of Keffer, Kaplan a. 
Yafet 

The elimination of R, and R, from system (13) leads always to the dispersion 
relation between w and X. Instead of using the secular equation, an indirect way can 
be taken, generalizing the complicated procedure applied in the paper of Keffer et 
al. (1953) (and also by Cofta 1956). 


From (13), putting for simplicity a = 0, P* = Q* = 0, it follows 


m: R,(f. A* + xB*) + Raf A” _ R, (AB? + xA”) + R,f,B* Ma 


Ry Re (37) 
_ (R,A" + RBI + (G,B* + RAD (« + fe) 
R, +R, 
Using an algebraic identity 
RA + RSB* = (R, + Ry) (4* + BY) — (RB* + R, 44) 
and introducing the abbreviation R,B* + R,A* = u, we obtain 
de dte La Ja Seine (38) 


R, +R, 


We can calculate the expression u also from (13), multiplying the first equation by R, 
and the second one by R, and afterwards comparing both sides: 


RiR, (4* — BY) (fı =)= (Ri B* — R347)f, 
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We multiply now both sides of this equality by (4* — B") and apply two identities: 
R,R,(A* — B*)? = (R, + Ry) (A* + B*) (R,B* + RA?) + 
— (RB + R,A*)? — (R, R) ABZ 
(RIB* — R2A*) (4* — B?) = (R, + R)? A*B* — (R,B* + RADA 


In this manner, we obtain for u a quadratic equation, equivalent to the secular equa- 
tion. The solution of this quadratic equation, substituted into (38), lead to the formula 


(14). 
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This paper presents the results of investigations of photosensitive counters used 
in direct reading spectrometry. ; 
Stable photosensitive counters constructed from metal with a cathode in the form of 

a strip set into the counter tube are described. The method used to measure the photo- 

electric efficiency is based on the known spectral distribution of an ultraviolet standard 

and the transmission of a UCI] 22 spectrograph. Counters with a platinum cathode 
proved to be the most stable. 
The variation of the photoelectric efficiency of G. M. counters was investigated. 

Three types of changes are discussed. y x 

a) Initial changes leading to a stable photoelectric threshold and a stable photoelectric 
efficiency. In order to abtain this about 10° counts are necessary. 

b) Reversible changes in efficiency of short duration are observed above a certain count- 
ing rate. These depend on the counting rate and the wave length, and occur with 
a certain delay. 

c) Changes in the photoelectric efficiency connected with the ageing of the counter. 
These changes involve a slow increase in the efficiency and a shift of the photoelectric 
threshold in the direction of the longer waves. They were absent only in counters with 
a platinum cathode. ! 

A relation between the slope of the characteristics for ultraviolet radiation and the 
photoelectric efficiency has been found: slopes between 5—55%/100 V correspond to 

efficiencies between 10-? and 10-19, 


Introduction 


One of the most sensitive detectors of ultraviolet radiation is a G. M. counter 
fitted with a quartz window. Although such counters have been known and used since 
1930 and many investigations have been devoted to them, there has been no systematic 
analysis of their properties. The reason for this is mainly due to the fact that the 
photoelectric effect in the presence of gases and vapours is very complicated and, 
in many details, not quite clear. The mechanism of the operation of a G. M. counter 
has not yet been sufficiently investigated as regards the influence of the photo- 
electric effect on the processes occurring during the discharge. 


(231) 
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One of the most important aspects of photosensitive counters which till now has 
not been satisfactorily explained is the question of the stability of the photoelectric 
efficiency. Karev and Rodionov (1935) tried to describe the connection between the 


photoelectric efficiency Y and the counting rate N by the relation Y= A + BYN, 
A and B being constants. These authors make no mention of the inertia of the changes 
in efficiency. 

By changing the length of the slit in the monochromator providing the light for 
the counters, Duffendack and Moris (1942) found that up to a counting rate of 
18000 cpm there is a linear dependance between the counting rate and the slit length. 
Katz (1954) also found that at low light intensities there is no relation between the 
photoelectric efficiency is independent of the counting rate for light quanta. In these 
three experiments non-selfquenching counters were used. 


Neuert (1948) has investigated counters using cathodes sputtered with silver 
and magnesium. He found that the photoelectric efficiency of a counter of this type 
when irradiated by a gamma ray source increases by a factor of the order of 10 and 
subsequently, in about 10 minutes after the gamma source is removed, the efficiency 
returns to its former value. This effect has been investigated by Neuert only for 
a wavelength of 2650 Å. | 


It is a well-known fact that counters with mica windows and glass counters begin 
to respond to daylight from incandescent lamps (A greater than 3000 A) after having 
counted about 10’ pulses. This matter is discussed in detail in the papers of Aron 
(1953) and Schwartz (1953). 


Neuert has measured large changes of the photosensitivity but he has not observed 
any changes for counters with brass cathodes. 

The present paper deals with cases in which the changes in efficiency were rather 
small. 

Another important problem connected with photosensitive counters is the ques- 
tion of the differences between the photosensitive characteristics for ultraviolet light 
and the characteristics for gamma radiation. Locher (1932) and Christoph (1936) 
pointed out that a greater photoelectric efficiency can be obtained by using a strong 
electric field at the cathode. 

This effect is caused by the back-diffusion of photoelectrons in the counter gas, 
which results in a considerable decrease in the probability that these electrons are 
caught by the electric field around the anode. The probability that a photoelectron 
will not return to the cathode increases with the electric field intensity E and decre- 
ases with the gas pressure in the counter p, i. e. it is proportional to V Elp. Christoph 
calculated that only about 10% of the photoelectrons do not return to the cathode. 
According to the most recent data (Theobald 1953 a and b) in pure argon, under 
counter conditions of pressure and electric field intensity, only 10 — 20 % of the photo- 
electrons do not return to the cathode. The probability of not coming back to the cat- 
hode increases with the decreasing energy of the photoelectron. At present there are no 
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other accurate data concerning back-diffusion, in particular no measurements have been 
made for monochromatic radiation and for the filling mixtures used in G. M. counters. 

Ito (1952), using a reduced hydrogen pressure in the counter, obtained a slight 
increase in the photoelectric efficiency. 

Wilkinson (1948) drew attention to the fact that the beginning of the plateau of 
the characteristics for ultraviolet is about 150 V above threshold, whereas for gamma 
radiation this voltage difference is only about 70 V. Wilkinson presented a theory of 
this effect based on considerations of statistical fluctuations in the initial avalanches, 
which can check the development of the discharge. The probability of stopping the 
development of the discharge along the wire decreases when the discharge in the 
counter is initiated not by a single electron as in the photoeffect, but by several elec- 
trons produced by nuclear radiation. 

The author performed systematic investigations on the problems described above 
and on the mechanism of counter operation. The aim of this work was, among others, 
the design of photosensitive counters for direct reading spectrometry. The results of 
studies on the mechanism of G. M. counters obtained by means of photosensitive 
counters are given in another paper (Ostrowski and Turek 1958). 


2. Apparatus 


The measurements were made on metal photosensitive counters with quartz 
windows, the construction of which has been described in several papers (Jurkiewicz 
et al. 1953, Ostrowski 1951, Buja et al. 1957). In Fig. 1 are shown two ways of mounting 
the quartz window by means of “Araldit“ 
manufactured by the Ciba firm of Swit- 
zerland. 

For the counter diameters used 


most frequently, outer dia. 16 mm, the 
maximum dimension of the windows 
were 10x 40 mm?. Usually windows 
measuring 5 X 20 mm? were used. The 
metal investigated had the shape of a 
bent springy strip. Its surface was cle- Fig. 1. Two methods of mounting quartz window 
aned with emery paper and washed with to the copper tube of a G. M. counter 
benzene. During soldering with a soft 

solder the counter was placed in a fixture in which running water was used to cool the 
window and the sample. About 80 counters were assembled, the majority having 
diameters of 16 mm. The counters were filled with an argon-alcohol mixture 
(120 mm Hg of argon and 12 mm Hg of ethyl alcohol). The counters used with gamma 
rays had a plateau length of 150 — 250 V, the slope of which was 1— 595/100 V. 
Counters with a platinum cathode had the smallest plateauslope. The counter back- 
ground was 0.8—1 cpm/cm? in agreement with the value given by various authors 
(0.5 —1.5 cpm/cm?). For counters with an Au and Pt cathode, which have a greater 
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gamma-ray efficiency, the background was 2 cpm/cm?. Counters with platinum 
cathodes had an efficiency for gamma rays from radium 2.5 times greater than copper 
cathode counters. 

The counters were tested at a counting rate of 200000 cpm at an overvoltage of 
100 V. After such tests it was found that counters with Zn, Cd, Al and aluminium 
cathodes became unstable (i. e. they flashed over). Al cathode counters withstood 
a test with a counting rate of about 100000 cpm, whereas Zn and Cd cathode counters 
already were burned out at a rate of about 50000 cpm. 

We can regard as stable, in the sense of this test, counters with cathodes of Ag, 
Au, C, Cu, Fe, Ni, Pt, brass, phosphor-bronze, copper-beryllium (295 of Be) and 
German new silver. 

The instability of the first group of counters occurs together with a shift of the 
photoelectric threshold towards the longer wavelengths for a period of the order of one 
hour. Thus, e. g. the threshold of Al cathode counters, after counting about 3 x 105 
pulses in a few minutes, shifts from 3200 À to the visible region. After the counter is 
disconnected, the photoelectric threshold shifts back towards the smaller wavelengths. 
The instability and the changes of the photoelectric threshold should be considered 
as closely related to one another. 

A special counter-box has been coustructed to be used in the UCII 22 quartz 
spectrograph instead of a plate box. The displacement of the counter together with 
the exit slit was made by means of a micrometer screw at a distance of 18 cm, so that 
the accessible range of the spectrum was 2030 —4500 À. A medium-pressure mercury 
lamp with liquid electrodes and an “Original Hanau S 300" lamp with optical proper- 
ties very similal to those of a standard ultraviolet source were used as light sources. 
The similar optical properties of these two lamps result from the nearly identical arc 
size in both. The lamp and the rate-meter used in the measurements took their power 
from a magnetic stabilizer. 


3. The measurement of the quantum photoelectric efficiency of G. M. counters 


The quantum photoelectric (or briefly, photoelectric) efficiency is defined as 
the ratio of the number of pulses produced in the counter by radiation of a given wave- 
length to the number of photons of this radiation falling on the cathode. Up till now 
thermopiles were used to determine the number of photons (Kreuchen 1935, Rodio- 
nov et al. 1955). This method has some disadvantages arising from the very large 
difference between the sensitivities of the counter and the thermopile. One is compel- 
led to use absorbers, reducing the intensity by several idees of magnitude, which can 
be the source of large error. 

Another method, based on calculations of the number of photons for a black 
body (Raich 1930, Katz 1954), has also been used. In this method one-can obtain in 
practice only very low counting rates, and the measurements can be anole, only for 
radiation of wavelength greater than 2600 À.. 

In the present work the known spectrum of the ultraviolet standard lamp was - 
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used and the transmission of the spectrograph in the various ranges of the spectrum 
have been calculated. The measurements of the spectral distribution of the standard 
ultraviolet source were made by Rössler (1939, 1952) by means of a quartz photo- 
electric cell calibrated with a thermopile. In this case the difference in the sensitivities is 
not so important as in the case of a thermopile and a counter. The standard ultraviolet 
lamp is a medium-pressure mercury lamp and it is known that the deviations of the 
relative intensities for particular spectral lines do not exceed a factor of 2, despite the 
different parameters of these lamps such as length of arc and quartz tube diameter (Kern 
1938, Sewing 1938, Meyer and Seitz 1947). Comparison with a standard source was 
of course made after reducing the intensity to the same input power (e. g. 100 W). 
From the work of the above-mentioned authors it follows that when using, instead of 
a standard source, a lamp of similar properties one obtains intensities for particular spec- 
tral lines which do not differ by more than a factor of 2, and probably considerably less. 
In the calculations of the spectrograph transmissivity one must take into conside- 
ration the radiation losses due to reflection from the quartz surfaces, the mirror 
surface of the aluminium collimator and also the absorption in the quartz. 
In order to neglect aperture losses the lamp should be placed at a suffi- 
ciently large distance from the entrance slit (in the present work this distance was 
1.4 m) along the axis of the instrument. The region in which the mercury arc. should 
be located can be determined experimentally by placing a strong lamp in the spectral 


transmision 
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Fig. 2. Transimission of the UCII 22 spectrograph. The data for the reflection coefficient of the alumini- 

nized mirror of the collimator were taken from the Critical Tables (1950, Vol. VI) and the data on the 

absorption of quartz, from Koller (1936) as cited in the book of Mayer and Seitz "Ultraviolette Strahlung" 
(1949). The scattering (about 596) was neglected 


plane of the spectrograph and observing a light spot (rectangular with cut corners) on 
a screen mounted in place of the staridard source. If we place the arc in this region, 
the radiation illuminating the slit will not illuminate any surface inside the spectro- 
graph besides the surface of the lens and the collimator mirror. By. using a wide 
enough slit we can avoid losses due to diffraction at the slit. The results of calculations 
of the spectrograph transmission are given in Fig. 2. 
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It is possible to de termine the transmission if we use two monochromators, so that 
the exit slit of the first is the entrance slit of the second. In this method we place the 
detector first behind the common slit and then behind the exit slit. Practically spea- 
king, this is a very difficult procedure because of the necessity of avoiding light losses 
due to the aperture. The results obtained by various authors (Katz 1954) for the two- 
monochromator method and for the transmission calculations are in agreement as 
regards the order of magnitude; differences of from 7 to 60 percent were noted 
in the various cases. 

It should be expected that the theoretical values should be greater than the 
actual transmission, owing to the scattering loss, which is taken as 5% for 10 layers 
(Angerer 1956 p. 261). 

The quantum efficiency of a photosensitive counter for the standard ultraviolet 
lamp at a distance of 1 m from the slit, is given by: 


EL k-On:S a) 


where N denotes the counting rate (in cpm) of photons corresponding to a given spec- 
tral line together with the background due to the continuous spectrum and natural 
background: 
Np, is the background due to the continuous spectrum and to the counters’ 
own background. 
k is the spectrograph transmission. 
O,, is the illumination in erg/cm? sec of the radiation corresponding to the 
particular ultraviolet standard spectral line at a distance of 1 m. 
S is the area of the entrance slit (in the present experiments 0.2 x 
x 0.0050 cm?). 
Thus O,,°S/hy is the number of photons entering the entrance slit in one second. 
For wavelengths below 2300 À no mercury lines were measured and we are obliged 
to use the continuous spectrum, the formula for which is similar to (1), namely: 
\Npa — Nal: hv 2 
ke Onesies (2) 


where Nz is the counters’ own background, 
Og, is equal to og, - 2: d,/sin a, 
where og, denotes the illumination in erg/cm?sec by radiation of wavelength A 
from the continuous spectrum of the standard source, at distance of 1 m, 
-z is the dispersion in A/mm, 
d, is the exit slit width in mm (0.2 mm), 
a is the angle between the axis and ta spectral plane. 

Since plane of the entrance slit is nearly perpendicular to the direction of the 
light (but not parallel to the spectral surface), it subtends a wavelength interval of 
z-d,/sina and k: Og, S is the power corresponding to the spectral interval sub- 
tended by the exit slit. 


Y= 
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If we use another mercury lamp at a distance of metres from the slit, we must 
multiply formulae (1) and (2) by M, - a?/M, where M, is the power of the standard 
lamp (250 W), M — the power of the mercury lamp which was enployed. 

From formulae (1) and (2) and from the results of Rössler the author prepared 
tables (for the line spectrum) and a graph (for the continuous spectrum) which per- 
mitted direct readings of the photoelectric efficiency for a given wavelength from the 
rate-meter readings. Over a period of several years the apparatus gave readings that 
were reproducible within the limits of 30%. 

The error was probably due to the mercury lamp with liquid electrodes, since the 
distribution of the mercury in the lamp was presumably random. 

The photoelectric efficiency obtained from the continuous spectrum was com- 
pared with that from the line spectrum. The both efficiency obtained from the 
platinum cathode counter was in good agreement. (The difference being smaller 
than 50%). Since the discrepancies between the results of various authors are rather 
great (several orders of magnitude) such an error should be regarded as sufficiently 
small. 

Recently, the results obtained with a mercury lamp S 300 (Original Hanau) 
were compared with those obtained earlier, and the results were found to be in agree- 
ment. As already noted, the properties of the S 300 lamp are very similar to those of 
a standard ultraviolet lamp. 

It should be remarked that the published comparisons of mercury lamps were 
made only for the spectral lines. The agreement between the efficiencies obtained for 
both the continuous and the line spectra-in the 2300 — 3000 A range allows one to 
suppose that below 2300 À one can also use the continuous spectrum to determine 
the photoelectric efficiency. Nevertheless, the deviation from the real values in the 
vicinity of 2050 A may be greater than that in the 2300 — 3000 A range. The cut-off 
in the measurements in the proximity of 2050 À is due to the end of the continuous 
spectrum of mercury. 


4. Results of the photoelectric efficiency measurements 


Fig. 3 shows the most typical variations of the photoelectric efficiency with 
wavelength for various metals. 


Counters with copper cathodes: 


The efficiency of the copper cathode counter shown in Fig. 3 was obtained after 
cleaning the cathode with emery paper. For the uncleaned cathodes the efficiency is 
about 30% of the value given in the figure. 

The reproducibility of the efficiency for various counter samples prepared 
under identical conditions is expressed by a factor of 2 for wavelengths below 2300 À. 
The use of various tubes and copper cathode elements of different origin did not 
E in any essential differences. The photoelectric threshold, taken as the maximum 
wavelength of the mercury line spectrum for which the counting rate is at least 100 cpm 


238 Kazimierz W. Ostrowski 


silver 


2000 2200 2400 2600 2800 3000 3200 
wave length 
Fig. 3. Results of the efficiency measurements of counters with different cathodes. A typical example 
was chosen for each sample. In the case of platinum the results for the same counter from May 1955 are 
denoted by the „+“ sign; from May 1956, by the „x “ sign. The measurements from May 1956 for the 
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Fig. 4. Comparison of the efficiencies of different counter samples with copper cathodes 
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in a given geometry was found to be within the limits of 2399 — 2965 A for copper 
cathode counters. The lack of a sufficient number of mercury lines did not allow this 
threshold to be determined more accurately. 

Fig. 4 illustrates the differences in efficiency for various counter designs. Counter 
No 5, was tightly sealed by "Araldit", and then heated at a temperature of 150°C 
for five hours in an atmosphere of air. 


Counters with brass and iron cathodes. 


Brass and iron cathode counters have a photoelectric efficiency distribution 
similar to that of the copper cathode counters, but the efficiencies are only about 20% 
of that measured for the copper cathodes. 

The reproducibility of the measurements was somewhat better than in the case of 
copper counters. 


Counters with cathodes from the alloy Cu + 2% Be -alloy. 


The interest in this alloy is connected with the fact that it is used as a material 
for electrodes in photomultipliers. The photoelectric efficiency of these counters was 
in some cases higher than that of copper cathode counters, even by a factor of 100 
(for 2800 A), but unfortunataly the reproducibility was not satisfactory especially : 
above 2400 A (e. g. by a factor of 10 or even more). The reason for this lack of 
reproducibility is probably connected with processes analogous to those, which are 
encountered in the activation of photomultipliers with Cu + Be -electrodes. This 
activation involves oxygen treatment at a high temperature. 


Counters with platinum cathodes. 


Counters with platinum cathodes and platinum iridium (2% Ir) cathodes have 
several very important properties in comparison with the other counters. These coun- 
ters had a higher sensitivity in the measurements described here at least in the spectral 
range below 3200 A. The reproducibility of the efficiencies in this case was rather 
good. Also the reproducibility of the threshold was the best for all types of counters 
investigated. Since the threshold is about 3500 A it can be determined rather accura- 
tely by investigating the radiation from an incandescent lamp, the slit of the spectro- 
graph and the counter being set wider. The slope of the characteristics of these counters 
for gamma rays is, as a rule, lower than 1%/100 V. In the first period after filling, 
changes in the Geiger threshold were observed to be greater than those in other 
counters. No satisfactory interpretation of this phenomenon was found. 


Counters with German silver cathodes (Cu + Ni + Zn). 


The photoelectric efficiency of these counters is only a few times smaller than 
that of platinum counters. This case shows the possibility of finding cheap alloys 
having good photoelectric efficiencies. The reproducibility of the efficiency was 
similar to that of the ^ platinum . cathode counters. 


240 Kazimierz W. Ostrowski 


Counters with cathodes made of silver, gold (22 carat), nickel and phos- 
phor-bronze. 

In certain wavelength intervals counters of these materials show an efficiency 
slightly greater than that of copper counters. The greatest efficiency was found for a gold 
cathode counter. Silver (Ag) and nickel coatings were obtained by electrolytical plat- 
ing. The symbol Ag’ in Fig. 3 refers to a cathode made from a strip of chemically 
pure silver. 

For other metals as aluminium, zinc, cadmium, lead, and tantalum no interesting 
results were obtained, i. e. their photoelectric efficiencies were in general much lower 
than those of the copper cathode counters, especially below 2400 À. 

Counters with cathodes of aluminium, zinc and cadmium were unstable. 


Carbon cathode counters. 


Carbon cathodes were obtained by evaporating India ink onto a copper strip: 
In the wavelength region above 2030 A the efficiency of this counter did not exceed 
2%, of a copper cathode counter. It was observed that in this case the sensitivity was 
associated with the copper itself and not with the carbon, i. e. a small part (2%) of the 
radiation is reflected by the area coated with India ink towards the uncoated copper, 
thus giving rise tothe effect mentioned above. 

India ink coating of the interior of counters has been successfully used in our 
laboratory to improve the stability of duraluminium beta-ray counters and of large 
aluminium cosmic ray counters. 


5. Comparison of the efficiency measurements with the results of other experiments 


The quantum photoelectric efficiency of photosensitive G. M. counters was 
determined in only a few cases. The measurements of Kreuchen (1935 a and b) and 
Rodionov et al. (1955) seem to be the most accurate. Kreuchen obtained the best 
results for zinc. This material, however, is not suitable for G. M. counters because of 
the instability of zinc cathode counters, especially at high counting rates. Almost the 
same photoelectric efficiency as that obtained by Kreuchen can be obtained for 
platinum cathode counters and with satisfactory stability. 

Rodionov et al. investigated counters with platinum cathodes coated by sputtering. 
These counters were filled with an argon-methyl mixture. The photoelectric efficiency 
of these counters below 2600 À is comparable to the efficiency of the platinum coun- 
ters described in the present paper, but for longer wavelengths it is lower (e. g. by 
2 orders of magnitude at 3100 À). The reproducibility of the threshold and of the 
efficiency given in the paper of Rodionov is smaller than that observed by the author 
(for about 35 cases). This might possibly be connected with the “maturing“ effect 
described in 6a. 

It is worthwhile to call attention to the paper of Katz (1954), who, in the interval 
2600 — 3200 A, obtained for non-self quenching counters with sputtered platinum 


n 2. 
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cathodes efficiences 30—50 times greater than those observed by the author. The coun- 
ters designed by Katz have a plateau of only about 50 V. Similarly, Labeyrie (1951) 
obtained an efficiency of the order of 10-? in the sensitivity maximum at 2500 À, 
for his counters with sputtered platinite (a Fe —Ni alloy with small additions of other 
metals). These counters were stable and had a good plateau. The efficiencies obtained 
by Katz and Labeyrie seem to be rather too high. 

The photoelectric efficiency of the metal, and especially the work function in 
vacuum after heating, have no simple connection with the corresponding quanties for 
counters containing gaseous mixtures. Gases soluble in metals have a strong influence 
on the photoelectric efficiency and on the work function. 

Platinum is a typical example; its photoelectric threshold in vacuum (2350 A) 
does not seem to encourage its use in photosensitive counters, while from the experi- 
ments described above, it follows that it is a suitable material for photo cathodes. 
This is probably connected with hydrogen ions being in solution with the cathode 
material. Another example is provided by the results of Eremin (1956) who obtained 
for polished magnesium and aluminium a threshold at 7300 A, where the counters 
showed a good photoelectric efficiency stability for several years, when a low counting 
rate was maintained. A change in the photoelectric threshold produced by the adsorp- 
tion of gases, is probably the reason for the different behaviour of some metals under 
conditions of high vacuum and under conditions met inside a counter. 


6. Three types of changes in photoelectric efficiency 


In the group of stable counters three fundamental types of efficiency changes may 
be distinguished: 


a) “Maturing” of of the counter. 

b) reversible changes of short duration connected with the process of counting. 
c) changes connected with the ageing of the counters during its operation. 
We shall now consider each of these changes in the photoelectric efficiency 


_a) “Maturing” of the counter. 

In order to stabilize the spectral distribution of the ne and. the. photo- 
electric threshold of a new counter, or one which has not been used'for a certain period, 
e. g., a month or more, it is necessary to count a certain number of pulses. For coun- 
ters with the internal diameter of 13 mm this stabilizing number of pulses is about 
(5 — 10 x 105 at an overvoltage of about 100 V. If the counter was “matured” pre- 
viously, then the required number of pulses is much smaller e. g. by one order of 
magnitude. Fig. 5 illustrates the results obtained for a platinum cathode counter. 
A definition of the photoelectric threshold entails basic difficulties, since the descrip- 
tion of the changes in the threshold are necessarily of approximate character. A certain 
indeterminacy of the photoelectric threshold occurring in vacuum is connected with the 
thermal motion. If the adsorbed substance has a definite potential barrier, then the 
indetaminacy of the photoelectric threshold will be further increased when the 
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tunnel effect is taken into account. The ,,maturing” of the counter in the sense men- 
tioned above is accompanied by a number of other changes (which can, at least partially, 
be reduced to photoelectric changes), such as the change in the slope of the characte- 
ristic, the increase in the plateau-length, the decrease in the length of the segment of 
the characteristic between the threshold and the beginning of the plateau. 


Efficiency 


2000 2200 2400 2600 2800 3000 3200 
wave length 


Fig. 5. Initial changes in the efficiency of a platinum cathode counter. After 6 x 105 counts at an over- 
voltage of 100 V no further changes were observed 
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Fig. 6. Plot of efficiency vs counting rate (in a state of equilibrium). The counter illustrated had a copper- 
berilium cathode (2% Be), internal diameter of 13 mm, overvoltage of 100 V. Measurement error about 3% 
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b) Reversible changes of short duration connected with the pro- 
cess of counting. 

If we measure the photoelectric efficiency of a counter at a low counting rate for 
a definite wavelength and if by bringing the gamma ray source closer to the counter, 
we increase the counting rate, e. g. to 50000 cpm for a few minutes, we observe, 
after removing the gamma-ray source, an increase in the photoelectric efficiency, 
e. g. by several percent. Changes greater than 5% have actually been measured. 

Investigations have disclosed some regularity in this effect. This is illustrated in 
Fig. 6. It was found that measurable changes of the photoelectric efficiency can be 
obtained only if a certain critical counting rate N, is exceeded. Let us denote the photo- 
electric efficiency of a counter below this critical value by Yo, and the efficiency at 
a given counting rate N by Y,(N). This efficiency automatically remains constant if 
the counting rate remains constant. The time sufficient for this is about 1 min. The 
efficiency measurement was made very rapidly (several seconds) at a low counting 
rate. In Fig. 6 is plotted the ratio Y,/Y, vs the counting rate N, for different wave- 
lengths. We see that the counting rate at which the above-mentioned increase in 
efficiency begins is independent of the wavelength of the radiation for which we are 
measuring the efficiency. Also, no relation was found between N, and the cathode 
material. No difference was observed for the two different filling mixtures (argon- 
alcohol, and argon-methylal) used in the work, although these last measurements 
obviously could not be very accurate. It was established that N, is at least approxi- 
mately proportional to the counter diameter and inversely proportional to the over- 
voltage. 

The following formula has been taken as an empirical working formula giving 
sufficient agreement with the observed data 


Y(N) = Y, E +A log a (3) 


The index r refers to the efficiency at "equilibrium". This formula is presumably not 
the only one which can describe the phenomenon. The introduction of the constant A 
is convenient, since it characterizes quantitatively the increase of the efficiency with 
the counting rate. 

Fig. 7 gives the slopes of Y,/Y 9 vs N curves for the values of the constants A 
in formula (3) for various wavelengths and materials used as counter cathodes. 

After removal of the radioactive source the photoelectric efficiency decreases. 
This decrease can be measured e. g. every 10 seconds for a period of one second, 
if we use the 0.5 sec time constant of the rate meter, by momentarily uncovering the 
window of the counter. It was found that Y — Y, decreases exponentially with time. 
'The time after which the efficiency Y drops to the arithmetic mean of the values Y 
and Y, is independent of the wavelength, the cathode material, and the counting rate. 
For argon-alcohol ensi it is 1.5 + 0.2 min and for argon-methylal mixtures, 


2.5 min. 
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For aluminium cathode counters, another time equal to (7 + 2) min has been 
found, in addition to the value 1.5 min mentioned before. In this last case the measu- 
rements were performed for 4350 A and 3650 A, i. e. for wavelengths for which the 
photoelectric efficiency is initially too small to be measured (< 107 12), In the case of 
unstable counters, such as those with aluminium cathodes, we can consider the whole 
phenomenon as a case with a variable Y, and a half-life of 7 min. 


brass 3x10 ? pulses 


Cu (metylal) 


new-silver 


2000 2500 3000 3500 
wave length 


Fig. 7. Variation of the constant A as defined by Eq. 3 with wavelength. As may be seen, the values of 4 
for the platinum cathode are the most advantageous. The least favourable are the values for an alumi- 
nium cathode. It should be noted that a small addition of Be to the copper is advantageous 


As already mentioned the decrease of the efficiency Y—Y, can be considered 
as exponential. Since the whole increase of the induced photoelectric efficiency does 
not last longer than 1.5 min there are experimental difficulties in investigating the 
nature of this increase. Let us consider as an example the data in Table 1, which 
gives the time in secs. required to reach the value Y, + Y,/2, starting from Y,. N, is 
the counting rate after attaining equilibrium. The counter used in this.test had a pla- 
tinum cathode. 


Table 1 


2399 T 
2965 14 
3120 18 
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As may be seen, the smaller the wavelength of the radiation used in the experi- 
ment and the higher the counting rate, the faster do we approach equilibrium. 

The greatest changes in the photoelectric efficiency can be determined in the 
following way. If N’ is the initial counting rate at which the photoelectric efficiency 
for some monochromatic illumination is Yo, the following relation expressing the pro- 
portionality of the efficiency and the counting rate holds: 

YY, = NIN? (4) 

where Y is the actual efficiency at the counting rate N (line b, Fig. 8). 


4 ry Y, 


7 2 3 4 5 6 7 8 9 
x 10 *pulses/min 


Fig. 8. Changes in photoelectric efficiency: a — equilibrium curve (Eq. 3), b, c, d — line, along which 
the change in efficiency takes place in the case of monochromatic illumination; A, B, C — sample 
initial points; D, E, F — final points of equilibrium. 


The intersection of line b with the equilibrium curve (equation 3, curve a), is the 
final equilibrium point (N,, Y,). Eliminating Y,/Y, from equations 3 and 4 by the 
method of successive approximations, we obtain N, as a function of N’/N,. 

If, in general, the photoelectric efficiency Y’ corresponds to the initial counting 
rate N’ i. e. if we proceed from an arbitrary point B on Fig. 8, then from the equation 


Y/Y, = NIN” (5) 


we can obtain a value N”, which corresponds to the efficiency Y, (line c, Fig. 8). 
Starting from the counting rate N” with the efficiency Yo, we pass through point B 
with the coordinates N’, Y’/Y,. This counting rate N” can be called the reduced 
counting rate. 

Fig. 9 shows the results of the N, calculations. The abscissa axis gives the devia- 
tions of the JV, — values from the reduced counting rate N” in %. The ordinates are 
the values of the ratio N”/Ny, i. e. N>- Yo/No: Y’. 
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Fig. 9. Changes in counting rate —Ne— as a function of ; » where N’ is the initial counting 
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Fig. 10. Changes in the photoelectric efficiency of a brass cathode counter as a function of the wavelength 
and age of the counter. The counter counted 100000 cpm at an overvoltage of 100 V 
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From this plot ye can determine the maximum error expected in using counters 
for spectral analysis, if the abscissae are greater than 1. The corrections for the dead 


time must be made separately. For the abscissa N’. Y,/N)- Y" < 1 (line d, Fig. 8) 
we obtain directly N, = N”, and therefore 


MU MM Ae M MEM: 
T Y! or AT — 2 Y? — py (6) 


The decrease in the efficiency in the region above the curve expressed by equa- 
tion 3 gives a different counting rate than the decrease along the line N = 0. The 
half-times for the decay of Y — Y are very difficult to measure more accurately. 
These values are somewhat greater than those in Table 1, but distinctly less than 
1.5 min. It is to be expected that if the decrease in the efficiency takes place at small 
values of N (in comparison with N,) we obtain times close to 1.5 min. 

The constant A increases with wavelength (Fig. 7). Platinum cathode counters 
have the most advantageous, i. e. the smallest, dependence of A on the wavelength. 
These counters have the lowest values of the A at the highest efficiency. It should 
be emphasized that the last point of the measurements for platinum cathode counters 
is only about 40 À from the photoelectric threshold (3500 À). The threshold in that 
case can be relatively easily determined by means of light from an incandescent lamp 
illuminating the spectrograph slit. 

A platinum cathode counter responds to daylight with a continuous increase of 
the counting rate up to the limit of its possibility (about 108 cpm). This effect is not 
an instability in the usual sence, because its response to gamma radiation or mono- 
chromatic radiation not too close to the threshold is quite normal. The response to 
continuous radiation in the vicinity of the threshold is in agreement with the relations 


given above (large values of A). | 
c) Changes in the efficiency and the constant 4 with the ageing 


of the counter. 

The dependance of the efficiency of a brass counter upon the life-time, expressed 
in terms of the total number of counted pulses is presented in Fig. 10. : 

For metals other than platinum the efficiency changes are very similar. The follo- 
wing characteristics have been observed: 

l. Changes in the photoelectric efficiency are not A E within the limits 
of the experimental error (about 20%) for wavelengths smaller than 2200 A. 

2. The photoelectric efficiency changes strongly during the first 2 x 10? counts 
and the photoelectric threshold shifts to a region where the glass and the mica are 
partly transparent (3000 A). This is the reason for the sensitivity of window counters 
to "light". 

3. The maximum photoelectric efficiency for radiation below 3100 À obtained 
in all cases investigated (Cu, Ag and brass) is approximately the same and is several 
times smaller than that found for platinum counters. The region below 3100 À could 
not be systematically investigated because of the lack of a sufficient number of spec- 
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tral lines of mercury and because of difficulties in the experiment for this region. 
It was established only that the photoelectric threshold drifts towards the visible 
region as the counter ages and that the photoelectric efficiency attained for this spec- 
tral region is small (of the order of 107%). 

4. The constant A defined by Eq (3), which in a measure of the increase in effi- 
ciency with the logarithm of thé counting rate, increases strongly with the age of the 
counter, e. g. by a factor of 10. No such changes were observed for counters with pla- 
tinum cathodes. 


7. Comparison of the results of the investigation of HEU stability with the 
results of other workers. 


The similarity between the Rodionov formula given in the introduction and for- 
mula (3) seems to be coincidental, since at the counting rates used by these authors 
there should not be any change in the photoelectric efficiency (Duffendack and 
Moris 1942, Katz 1954). 

The results of Neuert (1948) indicate large changes in efficiency. One may deter- 
mine from Neuert's data the half-life of the quantity Y — Y,. This time is about 2 mi- 
nutes. The value determined in this manner is in agreement with the value given in 
the present work (2.5 min for methylal). It seems that the counters with sputtered 
cathodes used by Neuert have a greater factor A than the counters used in the work 
of the author with massive cathodes or obtained by electrolytic methods (Ni, Ag). 
It would therefore be of value to compare sensitivity and stability of these two types 
of counters. The results of the work of Lauterjung (1947) on the regularity of the 
excitation and decay of efficiency concerns, we suppose, unstable counters, where 
these phenomena are still very complicated. The decay time for the difference of the 
order of 10 minutes is in agreement with the results of the present work for aluminium. 

One would expect to increase N, and decrease the factor A by using a fast quen- 
ching pulse (Picard and Rogoziński 1954). 


8. Counter characteristic for ultraviolet radiation 


The systematic investigations of the characteristics of photosensitive counters 
gave some new results on the existence óf a relation between the photoelectric sensi- 
tivity and the slope of the characteristic. 

a) Beginning of the plateau. 

Fig. 1l gives the characteristic of a photosensitive counter of 13 mm internal 
diameter with a platinum cathode. As may be seen, for gamma rays it is sufficient 
to raise the operating voltage of the counter by about 20 V above the threshold of 
the equipment (with sensitivity of 0.2 V) in order to attain the beginning of the plateau. 
For ultraviolet radiation the corresponding increase in voltage is about 50 —70 V. 
An unmatured, counter requires for gamma rays 2—3 times greater overvoltage in 
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comparison with a “matured” counter. This probably accounts for the difference 
between the results of Wilkinson and those of the author (see Introduction). 

b) Slope of the characteristic. 

The slope of the characteristic for a collimated gamma beam directed perpendi- 
cular to the counter axis is dependent on the point of incidence in the counter. In the 
case that the beam passes through the central part of the counter the slope of the 
characteristic is constant and is the smallest. This slope can differ from the slope of 
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Fig. 11. Characteristic of a platinum cathode counter for gamma rays and ultraviolet radiation. The basic ° 
difference occurs at the beginning of the characteristic 


the same counter irradiated by an uncollimated gamma beam by several percent, 
especially in the case of counters of small active length. 

The slope of the characteristic measured with a collimated gamma beam is a mea- 
sure of the number of spurious pulses, which depend on the voltage. The geometrical 
factor is eliminated in this way of irradiation. 

The measurement of the slope of the characteristic of G. M. counters for ultra- 
violet radiation entails some difficulties connected with changes in the photoelectric 
efficiency during operation. 

In order to eliminate the above-mentioned changes the counting rates were measu- 
red at two points of the characteristic alternately every 20 seconds. The measurements 
were made by means of rate meter with a time constant of 5 sec. A relatively low coun- 
ting rate was used — about 4000 cpm. It was found that the slope of the characteristic 
distinctly depends on the wavelength. By comparison of the slopes of the characteri- 
stics and the photoelectric efficiency of the counter the plot in Fig. 12 was prepared. 
It turns out, as seen in Fig. 12, that there exists a relation between the photoelectric 
efficiency and the slope of the characteristic of the counter for ultraviolet radiation. 
The slope for a collimated gamma beam passing through the centre of the counter was, - 
of course, taken from the slope measured for ultraviolet radiation. The deviation of 
the expetimental points from the curve does not exceed a factor of 5. The difficult 
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conditions of the measurement do not allow one, however, to expect better agreement. 
An error of 295/100 V in the determination of the slope of the characteristic is equi- 
valent to a factor of 2 in the photoelectric efficiency. For slopes of several ten percent 
(100 V an error up to 5%/100 V could be made. The error in determining the effi- 
ciency should not exceed a factor of 2. 

It follows from Fig. 12 that the greater the photoelectric efficiency the smaller 
the slope of the counter characteristic. A careful measurement of the slope for ultra- 
violet could even replace a measurement of the photoelectric efficiency. This method 
would allow a greater choice of radiation sources and a greater choice of the range 
of the radiation intensities used. 

The slope of the characteristics for ultraviolet radiations, in the opinion of the 
author, is connected with the phenomenon of the backward diffusion of photoelec- 


Efficiency 


10 20 50 40 50 %100V 
Plateau slope 


Fig. 12. Relation between the photoelectric efficiency and the slope of the characteristic for ultraviolet 
radiation @ and O represent counters with brass cathodes, other sings denote counters with a Cu + Be 
cathodes 


trons. This idea, in somewhat other form, was expressed by Locher (1932) and Chri- 
stoph (1936) (see Introduction). (The concept of slope of the characteristic was not 
yet known to these authors). __ 

Backward diffusion, in turn, is connected with the distribution of the velocity of 
the photoelectrons. The results shown in Fig. 12 represent, therefore, indirectly the 
relation between the photoelectric efficiency and the velocity distribution in the vici- 
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nity of the counter cathode. To explain this relation further detailed experimental 
data are necessary also in conditions of vacuum. Since the question of the relation 
between the photoelectric efficiency and the velocity distribution of photoelectrons 
is not discussed in the literature, this question should be regarded as an open one. 


9. Interpretation of the phenomenon of changes in efficiency 


The three types of changes in efficiency discussed above: *maturing", changes 
of short duration, and changes connected with the ageing of the counter are presu- 
mably due to the effect of the adsorbed molecules, mainly to their electric field. 

This effect may be of two types: 1) backward diffusion of photoelectrons in the 
adsorbed layer, and 2) tunnel effect and the phenomenon of quantum mechanical 
reflection of electrons in the electric field created by the adsorbed layer. 

Small changes in the shape of the potential barrier can evoke large changes in the 
probability of photoelectron emission, especially if the electron energy is smaller 
than the height of the potential barrier. In the vicinity of the photoelectric threshold 
the tunnel effect allows the determination of the probability of emission which is 
a known function of (E—U), where E is the energy of the photoelectron, and U the 
height of the potential barrier. If the energy of the photoelectrons is greater than the 
height of the potential barrier there takes place a quantum-mechanical reflection of the 
electron wave, which expresses itself in a coefficient of reflection. This coefficient is 
a known function of E/U. The continuous distribution of the photoelectron energy 
results in a certain fraction of the photoelectrons being transmitted through the bar- 
rier in a tunnel effect. 

A detailed calculation made by the author of the various heights and widths of 
the barrier gives changes in efficiency which are in agreement as to the order of magni- 
tude with that observed. 

Independently of the influence of the potential barrier cases are known in which 
an adsorbed substance of large dipole moment changes the work function of the 
metal in relation to vacuum. The change in the work function concerns phenomena at 
a depth of 10-8 cm and, in the absence of a potential barrier, has relatively small effect 
on the photoelectric efficiency far from the threshold. The structure of the deeper 
layers (at a depth of several tens of À) has the greatest effect on the efficiency. In the 
presence of the barrier, this influence is greater since there is still the above-mentioned 
probability of passing through the adsorned layer. For stable counters, it seems that 
the change in the work function does not occur over short periods of time. 

A strong dependence of the constant A on the wavelength which well accounts for 
the phenomena connected with the presence of a potential barrier suggests that the 
influence of backward diffusion in the adsorbed layer is slight in relation to the other 
effects; at least, such is the case for brief changes in the photoelectric efficiency result- 
ing from the counter operation. | 

The following picture may explain Eq. (3) given above. The threshold counting 
rate N, can be interpreted as a minimum one, or as a kind of "extrapolated" counting ` 
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rate necessary for the disassociation products to cover part of the cathode free from 
the adsorbed molecules of the quenching gas. Since these products, in the main, do 
not have a dipole moment, they do not give rise to a potential barrier. At a higher 
counting rate important changes take place, e. g. explusion of the adsorbed molecules 
of the quenching gas. The potential barrier then disappears in the regions in which 
this takes place. Other phenomena increasing the probability of photoelectric emis- 
sion are also possible. 

If the action of the new dissociation products were independent of the same 
products already adsorbed, then the relation 


YJY, = 1 +B (N—Ngl 


should hold. However this assumption is incorrect. The new dissociation products 
can displace the older ones. Thus the curve will be steep. 


10. Conclusion 


As a consequence of this work one may list some conditions which must be 
fulfilled in order for photosensitive counters to be used in spectrometry. 

1. The counters should be “matured”. 

2. Their operating voltage should exceed the threshold by at least 70 V. 

3. The choice of the cathode should be made not only with regard to the photo- 
electric efficiency, but also with regard to the changes in efficiency of short duration 
resulting from the ageing of the counter. There appears here a general characteristic 
of photosensitive counters, namely, a higher efficiency is accompanied, at least for 
a given material, by a better stability of the efficiency over both longer and shorter 
periods of time, and by a smaller slope of the characteristic. The behaviour of the 
counters with platinum cathodes is particularly favourable. 

' 4. The sensitivity and stability of the counters increase with the decrease in the 
wavelength of the light. The counter (with the platinum cathode) works particularly 
well below 2500 À, and sufficiently well below 3000 À. This behaviour is somewhat 
opposite to that of photomultipliers. 

It gives me great pleasure to thank Professor M. Miesowicz, Professor L. Jurkie- 
wicz, and Dr. J. M. Massalski for their critical comments and discussion. I would also 
like to thank the precision mechanics K. Biernacki and C. Czepielewski for their help 
in constructing the equipment. 
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SCATTERING OF SLOW NEUTRONS BY MOLECULES OF LIQUID 
WATER, AMMONIA AND HYDROGEN SULFIDE 


By A. WANIC 
2-nd Laboratory of the Institute for Nuclear Research, Cracow 


( Received December 18, 1958) 


Theoretical curves for neutron scattering cross-section vs. energy for H,O, NH; 
and H,S were compared with the experimental values for liquid substances. Disagree- 
ment obtained at low energy region was interpreted as an evidence of not existance of 
any free rotation of molecules in investigated liquids. 


1. Introduction 


The problem of neutron scattering by chemically bound nuclei was first treated 
theoretically by Fermi (1936). The main result of his considerations was the conclusion 
that proton cross-section exhibits the strongest dependence on the type of the bond. 

At present the best results in this field have been obtained by Krieger and Nel- 
kin (1957), who found the formula which is valid for all molecules and for neutron 
energy E, satysfying the condition 


E < Es = Eyib 


where E,,— is the energy of the rotational energy quantum of the molecule 
and E,, — energy necessary for excitation of the lowest vibrational level of the 
molecule. 

Krieger and Nelkin (K —N) have proved their theory in the case of gazeous 
H, — and CH, — molecules. The agreement with Melkonian’s (1949) experimental data 
was good within the above mentioned limits of neutron energy. Recent calculations 
performed by Janik, Maniawski and Rzany (1958) showed also an excellent agreement 
between the K —N theory and Melkonian's experimental data in the case of gazeous 
C,H, — molecules. 

Neutron scattering cross-section per proton for liquid water molecules, o£? (E), 
is experimentally known with a good accuracy since 1949 (Melkonian). Theoretical 
calculations of this cross-section were made by Danzker (1953) and Goriunov (1956), 
they did not however give agreement with experiment. Calculations based on K —/V 
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theory, performed by the author, also give an appearent disagreement (fig. 3). It 
is evident that the above discrepancies are due to the fact that none of the theories 
(including K —N) takes the molecular interactions into account. 

In view of these facts it was decided to measure neutron cross-sections of other 
simple liquids in order to compare the results with the theory of Krieger and Nelkin. 
As the subject of experiments liquid hydrogen sulfide and ammonia were chosen. 
Hydrogen sulfide's molecules are very similar to these of the water but are quite 
different in respect to molecular interactions. It could be suspected that liquid H,S 
fits better the assumptions of the K-N theory in the neutron energy range being taken 
into consideration. In ammonia intermolecular interactions show to a certain extent 
a similar character to those in water. 


2. Preparation of Samples 
Both, ammonia and hydrogen sulfide were prepared from reagents of high purity . 
using the reactions: 
NaOH + NH,OH > NaOH - H,O + NH, (1) 
NaS +H,PO, — Na4HPO,  --H,S(1) 


These reaction were performed in an apparatus shown in Fig. 1. In each of the cases 
mentioned, the liquid reagent was dropped, into the reaction vessel by manipulation of 


kKOH+CaO 


Fig. 1. Apparatus for preparation of liquid NH, 


the tap K. In order to avoid the contamination of gases by water vapor, drying sub- 
stances (Mg(CIO,), for HS and CaO-KOH mixture for NH;) and a cooled trap were 
applied. Cases were solidified in a hard-glass vessel, cooled by liquid air. 

As for neutron cross-section measurements it was necessary to have the liquids 
investigated in the form of a layer of well defined thickness, special metal vessels, made of 
aluminum-alloy, were constructed (Fig. 2). The liquid containing chamber was made 
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by screwing together two parts A and B, smearing the screw-threads with araldit- 
-glue. After polimerization of araldit the vessels were examined under 60 at. pressure. 
Then they were filled with liquids (H,S, NH,) at an appropriate low temperature, and 


Fig. 2. Aluminum vessel for transmission measurements 


closed with steel cones pressed down into holes h, and h,. At room temperature the 
pressure in the vessels was about 10 at and 20 at for NH, and H,S respectively. 

On the next day the vessels were photographed using X-rays. On skiagrams the 
situation of meniscus could be seen. 


3. Cross-Section Measurements 


` Proper measurements were perfermed using reactor neutrons, but some intro- 
ductory investigations for H,S were also made with Maxwellian neutrons. The source 
of neutrons consisted of 100 mg of Radium mixed with Berylium powder placed in 
the centre of a paraffin wax block. With such a source it was not possible to work in 
a good geometry. This was the reason that not absolute measurements but only relative 


to water were performed. 
Ratio of H,S to H,O cross-section was calculated from the formula 


OHS . HO My,s %H,0 A Mny,o In Iu,o (1) 
Oy,0  Cns/muoXms  0H,S %H.S H,O Ins 
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where 9 — is the density of the investigated sample, 

x — thickness of the sample, 

M — mass of molecule, 

I — the intensity of neutrons passing through a layer x 
and ø — total neutron cross-section 


Ino 


x was chosen so that In 
H,O Ius 


<1 (i.e. I o ~ Ig, s). Water was placed in a similar 


aluminum vessel as that for H,S. Transmisions of the empty vessels were compared in 
this same geometry in which measurements were made. The transmisions were iden- 
tical within the limits of statistical error. 

Measurement of neutron cross-section gave the result: 


PHS _ 0.963 + 0.015 


94,0 
This value was ascribed to neutrons having energy equal ca 0.04 eV (1949, 1941) in 
satisfactory agreement with later results obtained with monochromatic neutrons. 


The main experimental work was done with the monochromatic neutron beam, 
after the Polish reactor “Ewa” was put into operation. A simple crystal monochro- 
mator was mounted near the output of one of its channels. It consisted of an aluminum 
monocrystal and an arrangement to turn the crystal around the vertical axis. The energy 
of neutrons reflected by the crystal was not calculated from the Bragg equation but 
was obtained from measurements of transmission in a silver absorber taking into 
account the Ag-cross-section tabulated data. 

Due to a satisfactory intensity of the neutron beam it was possible to measure 
the cross-sections in an usual absolute manner using the formula: 


M 
FE YE (2) 


I : 
where T — 72 denotes the transmission of neutrons in the investigated liquid. 
0 


I, was measured when the beam passed through the empty aluminum vessel. Densi- 
ties of liquids were obtained by interpolation of tabulated values to room temperatures. 

Cross-sections of H,S and HN, were measured for seven different neutron ener- 
gies. In the measurements with three lowest neutron energies a beam-filter (thin 
cadmium foil) was used in order to decrease the beam contamination in neutrons of 
second order reflection. 

The background was detected by cutting off the direct beam with a thick cad- 
mium foil (1 mm). 

In all investigations as a detector of neutrons a BF, proportional counter was 
used. i 

The number of counts required for a necessary statistical accuracy was gathered 
in repeated single measurements with succesive interchanging of samples. This process 
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had to cancel the errors which could have been introduced by slight instability of the 
reactor power level or of the parameters of electronic equipment. Experimental results 
are shown in the fig. 4 and 5. 


4. Theory 


Theoretical calculations of the ~.oss-sections were performed using Krieger and 
Nelkin theory. According to that theory scattering cross-section of proton bound in 
a free molecule is given by the formula: 


ori (En) = Stree g (erf [ot] — (1— pt e7 erf [eè (1 — p)} 8) 


al! m --M 2M 

»--5 NE i iim 

T is the temperature of the semple, 

M — effective mass of scettering proton, 
and œ — its effective frequency of vibration. 
M was calculated taking into account the following data: for H,O: A, and A, — values 
given in Krieger and Nelkin's paper, for H,S: molecular structure constants measured 
by Burnus and Gordy (1953) 
and for NH,: structure constants quoted by dbz (1956). 

In order to obtain w, amplitudes of normal vibrational modes in a given mole- 
cule must be known. In the present calculations they were got from a known geometry 
of molecular vibrations (9) and a normalization condition. By aubstituting w ‘and 
M — values calculated for each molecule to (3) following results (for kT = 0.0253 eV) 
were obtained: 


13.93 - > 
oio (En) = — p (erf (9.0086 YE,) — 0.9841 e -2:56 En erf (8.865 VE,)} (4) 


10.21 * ! 
ons (Ex) = —p— (erf (9.485 V E) — 0.9792 e = 3.697 En erf (9.287 LEE (5) 


7.880 
oN (E) = ins (erf (10.613 |/E,) — 0.9758 e- 5-394 E, erf (10.356 VE,)} (6) 


In addition, some NAME calculations were performed. 

In the case of H,S it was assumed that the sulphur nucleus has an infinite mass 
possesing however unchanged parameters of molecular structure and unchanged 
frequencies; assuming that one gets: 


E 2l 


(> (Ej)]e = us (9.6904 VE.) — 0.9797 e- 3-766 E, erf (9.494 VE,)} (7) 
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In the case of water, contemporary knowledge of the structure of its liquid setat 
was taken into account. There are facts from which the conclusion may be drawn 
that the molecules in water do not rotate freely but are torsionally oscillating (10). 
The frequency of this oscillation is about 500 cm. Hindrance of rotational motions 
was introduced in equation (3) by taking M= Mno = 18 a. m. u., and v* was treated 


formally like the other molecular frequencies, reducing w. In this way was obtained 


zx 
[02° (E)]rauia = 


eu (26.54 V E) — 0. 9826 e— 24.27 E, erf (26.08 VE) (8) 


Of course the validity of the above formula must not be expected for E, > hv* ~ 
Results of the calculations described above are shown on the Figs. 3,4 and 5 
together with the experimental data. 


5. Discussion of Results 


The results presented in preceeding sections show that for the liquids investigated 
molecular cross-sections are greater than the calculated for free molecules ones. This inc- 
rease of cross-section (connected with the liquid state of the scatterer) may be interpre- 
ted as an evidence of complete or partial hindrance of free molecular rotation in liquids. 


Oy 
(barns) 
55 


50 


2 an i UM 0.08 010 012 En(eV) 


Fig. 3. Neutron cross-section per proton in H,O-molecule vs. energy 
= one (E,) for free molecule; theory, 
—— ee (E,)] liquid; theory, authors calculations 
a oe [cto (E,)] liquid; theory, McReynolds and al. DOSE) 
O Melkonian's experimental points for liquid H,O. 
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0.02 0.04 0.06 0.08 070 0.12 Enfev) 
Fig. 4. Neutron cross-section per proton in H,S-molecule vs. energy; 

— aHıS (E,) for free molecule, theory; 

— — — gas (E,)theory with the assumption: M, = 2o 

© experimential points for liquid H,S. 


0.06 0.08 0.10 012 En(ev) 


002 0.04 


Fig. 5. Neutron cross-section per proton in NH;-molecule vs. energy 
— (o y Pa) for free molecule, theory; 
© experimental values for liquid NH. 
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The hindrance of free rotation in water seems now to be rather well established 
by optical and neutron experiments (Brockhouse 1957, Cartwright 1936). Recently 
this fact was confirmed by the work of McReynolds et al. (1958). They proved that 
the behaviour of neutron crosssection of water in a low energy region is explainable 
by a liquid-water model with a hindered rotation with Àv* = 0.06 eV. Water cross- 
-sections evalusted by then are represented in figure 3. A divergence from the 
author’s results is connected with some differences in approximations. 

Data on intermolecular forces in ammonia (evidence of strong hydrogen bond) 
imply that mechanism of molecular motions in this liquid is more or less like this in 
water. 

A generalization of this picture for liquid hydrogen sulfide is however unjustified 
by the present state of knowledge. A model of liquid H,S in which translational mo- 
tions are completely frozen only slightly corrects the results of calculations (Fig. 4). 
The hindrance of rotation, which must be therefore accepted also in this case, seems 
not to be caused by specific features of atractive forces (hydrogen sulfide does not 
show a tendency to association by hydrogen bonds) but is probably connected with 
geometrical assymetry of closely packed molecules. 
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NOTE SUPPLÉMENTAIRE A L'ARTICLE: 


“L'influence des moments quadrupolaires sur la constante diélectrique“. 
A. R. FERCHMIN 
Institut de Physique Théorique, Université Adam Mickiewicz, Poznan 


Dr A. D. Buckingham a bien voulu attirer l'attention de l'auteur sur la 
nécessité de considérer l'action polarisatrice du champ de réaction produit par le 
moment quadrupolaire sur le dipole. En effet nous obtenons ainsi dans la formule 
pour la polarisation un terme dépendant des moments. dipolaires et quadrupolai- 
res comparable au terme calculé par l'auteur (A. R. Ferchmin, Acta Phys. 
Polon., 18, 133 (1959). Il jouera quelque róle dans le calcul des corrections 
proposées par l'auteur dans le travail mentionné ci-dessus, mais ne changera pas 
les conclusions qualitatives concernant le signe du moment quadrupolaire de la 
molécule. ; 
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LABORATORY EQUIPMENT AND TECHNIQUES 
A UNIVERSAL DOUBLE-CRYSTAL NEUTRON SPECTROMETER 
By Denis O’CONNOR AND LucjAN BoNKOWSKI 


Institute of Nuclear Research, Polish Academy of Sciences, Warsaw 


(Received October 23, 1958) 


The construction of a universal double crystal neutron spectrometer of high preci- 
sion is described. The instrument can be used for double or single crystal neutron spectro- 
metry, as a neutron monochromator for cross section measurements in the range 0.03 to 
10 eV, and as a diffractograph for powder diffraction crystallographic studies. Some 
illustrative results of measurements are given. 


The crystal spectrometer to be described below was designed for use in the follow- 

ing applications: (a) as a double or single crystal neutron spectrometer for neutron 
spectrum studies and neutron cross section measurements in the range 0.03 eV to 
10 eV, (b) as a diffractograph for high precision crystallographic studies. Several 
descriptions of instruments of this type are to be found in the literature (1). 
. The design of the instrument was required to fulfill the following conditions: 
(i) large neutron beam cross section and high transsmision efficiency (ii) universal 
in application i. e. provision of two crystal tables on the same instrument to enable 
simple rearrangement for applications (a) and (b) above, (iii) electrical drive and 
remote control of principal movements, (iv) high accuracy (1) of measurement and 
high reproducibility of angular positions of crystal tables and main spectrometer 
arm. To fulfill condition (i) above requires an instrument of considerable size to acco- 
modate the large monochromating crystals, collimators and counters to give the 
required transmission. Condition (ii) implies the use of a second auxiliary crystal 
table, mounted off the main spectrometer arm and table axis, and simple means for 
moving the whole base of the machine perpendicular to the neutron beam to bring 
either of the crystal tables axes into the beam. Condition (iv) is fulfilled by suitable 
mechanical design and to some extent by accurate machining of some of the compo- 
nents. 

The main features of the instrument will be appreciated from Figure 1, 2 and 
the photograph. The important dimensions and tolerances are presented in the table. 
The most important mechanical feature, which to a large extent determines the rest 
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of the design, is the method of support of the main spectrometer arm. This arm, 
whichis 2.5 metres long and will carry an evenly distributed load of 200 kg, is supported 
at about its centre of gravity by two ground steel wheels running on a steel rail formed 
on two curved beams forming part of the base of the instrument. Using this construc- 
tion the height of the arm above the floor has been reduced, the weight of the machine 
held reasonably small (no counter weights for the arm are neeed) and the plane of 
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Fig. l. Plan diagram of spectrometer. 


K, — collimator in experimental hole 

R, — main spectrometer arm 

R, — auxiliary arm carrying crystal table Ss 

Sı — central crystal table 

S, — sliding supports on arm for collimator and counter; 
S, — automatic sample holder 


the arm maintained with good accuracy — under constant load the deflection of the 
tip of the arm varies only by 0.4 mm in 120° rotation. The arm is rotated by driving 
one of the steel supporting wheels. The elastic distortion of the rail by the wheels is 
reponsible for reducing the overrun of the arm to such an extent that it stops within 2 
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seconds of arc when the drive is switched off. The upper surface of the arm is provided 
with machined ways along which slide machine tool-type supports for the counter 
and collimators. This facilitates simple adjustement of the position of the latter. 

The curved beams and rails form part of the three — legged base of the instru- 
ment which also runs on steel rails set in the concrete floor. The base can be levelled 
to bring the spectrometer arm in the horizontal plane with an accuracy of 0.02 mm 
in 1000 mm and this level is maintained with movement of the instrument across the 
neutron beam. 

The central crystal table which is mounted on the main arm is coupled to the 
arm by a steel pulley and tape system to give a ratio of table to arm speed of 1:2. 


This ratio is maintained to 0.1’ in 60° arm rotation. Means is provided for rotation 
of the table independently of the arm, and without breaking the coupling. The surface 
of the crystal table is 160 mm below the optic axis, this providing ample room for 
a variety of goniometers and crystal holders. 

The auxiliary crystal table is mounted on a short arm (60 cms between table 
axes), rotating in bearings fixed to the main spectrometer arm, this holding the axes 
of the tables parallel idependently of the deflection of the main arm. This ensures 
that the counter sees only those neutrons which pass through the axes of both crystals. 
The second crystal table is finished off in an identical CES to the first so that the 
crystal holders fit it. i 

The angular poa of the main arm, the first NS table, and the auxiliary 
arm are measured to l'with an accuracy of 0.1’ by means of divided glass scales and 
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optical projection systems. Simultaneous measurement of the angular positions of 
the crystal table and the main arm provides a check of the table arm coup ling ratio. 
A selsyn fixed to one of the wheels of the spectrometer arm provides an additional 
means of remote angle measurement with an accuracy of l'. 

The Soller collimators consist of steel strips clamped together in a steel block, 
with regular spacing between them to form a number of long thin channels of the 
required angular width. The cross section of the emergent beam is 4 x 4 cm2. Both 
collimators k, and k, are of identical construction, k being made to fit the experimental 
hole. Two sets of collimators are in use of half angle 10’ and 20’. 

The counter used is a boron trifluoride proportional counter made in the Electro- 
nics Laboratory of this Institute, of 47 mm diameter, 720 mm effective length, and 
filled to a pressure of 600 mm Hg. The end window of the counter is of copper 1 mm 
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Fig. 3. Rocking curve for copper single crystal. Reflection of 1A neutrons off 111 plane. Collimator 
angular width 10' 


thick. An enriched boron counter is used to improve efficiency for higher neutron 
energies. The counters have a plateau of 300 volts, of slope 295/100 V. A second small 
counter placed near the exit of the first collimator is used to monitor the intensity of 
the beam. 

All the electrically driven movements of the spectrometer can be controlled 
remotely from a single switch panel. This panel together with the pulse amplifiers, 
scalers etc. are mounted at a distance of about 10 metres from the spectrometer. 

The measurement cycle of the apparatus is completely automatically controlled, 
and the results of measurements automatically recorded by apparatus which has 
‚been described elsewhere (2). 
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A typical example of the measurements made on this machine is illustrated in 
Figure 3. This shows a rocking curve for a copper monocyrstal (3) for the 111 plane 
reflection for 1 A neutrons. The collimators used in this case were of 10’ half width. 
The very small increase in half width due to the crystal indicates that the latter has 
a very small mosaic spread. At a reactor power of 200 KW the counting rate of the 
counter is about 140.000 counts per minute for this crystal at 1 A neutron wavelength. 
At this wavelength with 10’ collimators the wavelength resolution is of the order 
of one percent.. 

The authors are indebted to Prof. Bronislaw Buras for his advice and encourage- 
ment throughout the course of this work, to ing. W. Ney and the staff of the Design 
Office of the Institute of Nuclear Research for advice and cooperation in the design 
stage, and to the technical staff or the "Rafamet" Machine Fool factory, Racibórz, 
for their care and skill in the construction of the spectrometer. 
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Table 


Characteristic data of neutron spectrometer 
Main spectrometer arm: 


length 2500 mm 
angular speed with electric drive 5 or 40 per min. 
power consuption of drive motor 12 W 
accuracy of angular measurement 0.1 
maximum loading 200 kg 
angle of rotation —60° + 110? 
Crystal tables: 

diameter 200 mm 
table — optic axis distance 160 mm 
angular speed with electric drive 5 per min. 
accuracy of angular measurement 0.1 
Accuracy of positoning of main axis 0.5 m 


Total weight 1500 kg 
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